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Abstract 


The impact of land-surface properties like vegetation, soil type and soil moisture, 
and the orography on the atmosphere is manifold. On one hand, these features 
determine the evolution of the atmospheric boundary layer, in particular, the con- 
vective conditions and further the exchange of mass, momentum, heat, and hu- 
midity with the free troposphere. Subsequently, the land surface also influences 
the pre-convective environment. On the other hand, land-surface heterogeneity 
results in the spatial variability of the land-surface parameters, which often leads 
to thermally-induced circulations and associated convergence zones in the con- 
vective boundary layer. These, in turn, act as trigger mechanisms for convective 
clouds and precipitation. That means in simulations, the distribution and amount 
of clouds and subsequent precipitation depend on the resolution of land-surface 
and model grid spacing (A,) alike. Therefore, the focus of the study is to (1) com- 
pare areal mean precipitation for different model grid spacings and land-surface 
resolutions, (11) analyse reasons for their differences, (111) investigate spatial pre- 
cipitation patterns, and describe relevant trigger mechanisms of convection. 

The impact of model grid spacing and land-surface resolution on convective pre- 
cipitation is investigated within the framework of the HD(CP)* project (High 
Definition Clouds and Precipitation for advancing Climate Prediction). For that 
purpose, geographical areas with different types of complexity in the orography 
and considerable number density of lightning strikes (deep convection) are selec- 
ted. The areas are: the flat terrain near Berlin (A1), the isolated Harz mountain 
range in central Germany (A2), and the complex terrain, the Black Forest moun- 
tains (A3). Six suitable days with weak large-scale forcing but a considerable 
number of lightning strikes are chosen. ICOsahedral Nonhydrostatic (ICON) sim- 
ulations in large eddy model setup have been performed using six model grid 
spacings: Numerical Weather Prediction (NWP) mode (Asooom; A2500m), Large 
Eddy Simulation (LES) mode (Aj259m, A625m, A312m and Aj56m) in a nested do- 
main setup (control runs). The 41567 control run is the reference run. The impact 
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of land-surface resolution on areal mean precipitation and precipitation patterns 
has been deduced by reducing the resolution of land-surface properties, e.g. ve- 
getation, soil type, and the orography. The differences in simulated areal mean 
precipitation are explained through heat and moisture budget calculations. Vari- 
ations in the precipitation patterns are analysed by investigating relevant triggering 
mechanisms. The source regions of the convective precipitation are identified by 
applying a backward trajectory model. To diagnose the turbulent sensible and 
latent heat fluxes at the Earth’s surface in the source regions of convection, their 
dependence on parameters like orography, soil moisture index, transpiration area 
index, and net radiation is determined using the standardised multiple regression 
techniques. 

The results show that the areal mean accumulated precipitation amount for 
most of the cases decreases systematically across the LES grid spacings from 
A1250m to 61567. The relative precipitation difference normalised by the precipit- 
ation in the reference run is in the range of -26 to 400 % with the 75th percentile of 
155 %. In four out of the six days, 6125072 results in intenser precipitation patterns 
and an earlier onset of precipitation by | to 2 hours in comparison to the reference 
run. The modification of land-surface resolution from 156 m to 1250 m leads to 
variability in the mean precipitation in the range of 17 to 37 % with the 75th per- 
centile of 7 % which increases to a range of -17 to 49 % and the 75th percentile 
of 22 % with the land-surface resolution of 5000 m. The land-surface sensitivity 
experiments show a negligible impact on the onset time of precipitation and the 
precipitation patterns. Thus, the modification in land-surface resolution results in 
much smaller variability in the areal mean precipitation amount in comparison to 
the model grid spacing. 

To understand the differences of areal mean accumulated precipitation and 
onset of precipitation between the control runs, the heat and moisture budgets are 
analysed in detail for one day, for which the relative difference in the mean precip- 
itation by Aj259m and 61567 is ~175 %. Unlike Aj259m, A156m first shows intensive 
evaporative cooling due to the formation of numerous small clouds. Evaporative 
cooling is generated at the edge and shell regions of the small clouds. As a res- 
ult, the clouds often dissolve before they could grow deep enough to precipitate. 
In the subsequent hours cloud aggregation is a crucial step causing precipitation 
generation in 61567. 
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Concerning initiation of convection, overall the LES grid spacings show the sim- 
ilar thermally- and orographically-induced circulations in all areas (Al, A2, and 
A3). However, as demonstrated for Al considerable differences in triggering 
could occur when the land-surface resolution is reduced down from 156 m to 
5000 m. This finding holds when the resolution of 5000 m smoothed out those 
land-surface heterogeneities (e.g. lake breezes and urban heat island) which are 
responsible for convection initiation at 156 m land-surface resolution. 
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1. Introduction and motivation 


Clouds are complex bodies because they consist of numerous droplets interact- 
ing through a variety of cloud-physical processes. For example, its growth and 
dissolution due to evaporation are very complex and take place at a very fine 
scale. On the global scale, clouds are one of the largest uncertainties in the 
climate projection by models (Bony and Dufresne, 2005; Boucher et al., 2013) 
and NWP models (Zhang et al., 2019). Numerical models provide opportunities 
to study the different multi-scale processes and their interactions. As the NWP 
models at high resolution can resolve the deep convection, the parametrization of 
deep convection is not required for model grid spacings of a few kilometres and 
finer (Weisman et al., 1997). The convection-resolving simulations of short-range 
NWP (Weusthoff et al., 2011) and longer-range projections (Hohenegger et al., 
2008), explicitly resolve the governing nonhydrostatic processes of deep convect- 
ive storms. Hence, they perform better than convection-parametrizing simula- 
tions. However, the resolved convection and the resulting precipitation on kilo- 
and hectometre scale may vary depending upon the embedded resolved (and un- 
resolved) processes. 

Therefore, the selection of suitable model grid spacings to study deep con- 
vection, clouds, and precipitation is a difficult task as it poses questions about the 
extensive usage of the computational power and the complexity of resolved and 
other parameterized processes like turbulence and cloud-microphysics. Convect- 
ive precipitation is one of the difficult phenomena to capture in the model sim- 
ulation. There are several challenges in the forecast of precipitating convective 
events, e.g. the onset time and the precipitation patterns that could be dependent 
on the model grid spacing, parameterization schemes, initial and boundary condi- 
tions, land-surface resolution, etc. Model simulations are sensitive to the factors 
mentioned above; therefore, it should be chosen carefully based on the applica- 
tions. With this, the forecasting systems using kilometre to hectometre resolution 
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provide insights into the scalar dispersion, city ventilation, cloud modelling and 
wind energy production (Liu et al., 2011; Zhu et al., 2010; Taylor et al., 2016). 
Concerning moist convection, the processes of Convective Boundary Layer (CBL) 
are quite decisive. The CBL properties determine the convection initiation, which 
are, e.g. represented by parameters like Convective Available Potential Energy 
(CAPE) and Convective Inhibition (CIN). The spatial heterogeneities of the CBL 
often generate mesoscale circulations which then act as trigger mechanisms for 
convection. The state of the CBL, in turn, is determined by the surface fluxes, 
1.e., how the available energy at the Earth's surface 1s partitioned into the sensible 
and latent heat flux. This partitioning in turn depends on the land-surface condi- 
tions (Kalthoff et al., 1999; Western et al., 2002; Koster et al., 2004). Land-surface 
heterogeneity results from the variations in land cover and land use across a range 
of length scales, 1.6. the evolution of moist convection are strongly coupled with 
the characteristics of land-surface types (Avissar and Chen, 1993). 

The aforementioned factors and processes influence the distribution and amount 
of convection, clouds, and precipitation. But the question of characteristic model 
grid spacing in conjunction with an optimal land-surface resolution to study the 
moist convection processes is still difficult to address as their relative sensitivity 
to convective precipitation over scale of a kilometre and hectometre is not well 
known. With this overview, the scientific objectives addressed in this thesis are as 
follows: 


1. To investigate the impact of model grid spacing and land-surface resolution 
on clouds and precipitation. 


2. To monitor the processes causing different precipitation behaviour across the 
model grid spacings. 


3. To identify the triggering mechanisms over the grid spacings resulting in 
different precipitation patterns. 


To achieve the above goals the ICON model has been used in this thesis which 
is a crucial motivation of this thesis and allows to test the capability of the new 
ICON model. This study focuses on the inter-comparison of locally induced con- 
vective precipitation simulated in ICON over a range of model grid spacings and 
land-surface resolutions. ICON is a unified model system which is now oper- 
ational at Deutscher Wetterdienst (DWD) since January 2015. This model has 
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been developed in collaboration with DWD and Max-Planck-Institute for Met- 
eorology (MPI-M). It can be run in three different modes depending upon the 
application of the users. These modes are: ICON-ECHAM, ICON-NWP, and 
ICON-Large Eddy Model (LEM). The details of the general configuration and 
model setup are discussed in Chapter 4. ICON atmospheric model in LEM mode 
has been used in this study. LEM is a setup in limited area mode with LES turbu- 
lence. 

Convective precipitation is often localised, intense, and short-lived in nature. To 
identify the hot spot areas of deep convection, the number density of lightning 
strikes for the summer period (MJJAS) from 2000 to 2018 is investigated, light- 
ning being a proxy of deep convection (Leary and Ritchie, 2009). The lightning 
data is retrieved from Siemens lightning information service (BLIDS), which is a 
part of the European Cooperation for Lightning Detection (EUCLID). The details 
of the used lightning data are discussed in Chapter 3. Deep convection preferably 
occurs over orographically complex terrain (Kottmeier et al., 2008; Kalthoff et al., 
2011; Barthlott and Kalthoff, 2011). Therefore, depending upon the complexity of 
the orography and the associated number density of lightning strikes, three areas 
turn out to be most suitable which are: the flat terrain near Berlin (Al), the isolated 
mountain range, Harz mountains (A2) and the complex terrain, the Black Forest 
(A3). 

As the impact of land-surface resolution is expected to be most prominent for 
locally induced circulation, the conditions with weak large-scale forcing over the 
above-selected areas are chosen. In total six cases are chosen, two for each area, 
in this study (Chapter 3). These days have been analysed by performing the ICON 
simulations (See Chapter 5 for details of the configuration of ICON simulations). 
An overview of the relative impact of model grid spacings and land-surface resol- 
ution on convective precipitation is presented (Chapter 6). The reasons for differ- 
ences have been analysed using the heat and moisture budgets. The budget com- 
ponents have been implemented in ICON in online mode. The budgets infer the 
role of heat and moisture in regulating the processes of mass, momentum, heat and 
moisture exchange between the Earth's surface and CBL, and further to free tro- 
posphere activated at a particular grid spacing. Furthermore, to identify the source 
areas of this triggered moist convection, the backward trajectories are calculated 
using the trajectory model LAGrangian ANalysis TOol (LAGRANTO) (Wernli 
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and Davies, 1997). LAGRANTO has been adapted for ICON model fields (for 
details of LAGRANTO see Section 5.6). The backward trajectories facilitate the 
understanding of the conditions of the source region resulting in the formation of 
convective clouds (Chapter 7). 

To understand the impact of land-surface resolution on precipitation, the sens- 
itivity simulations have been performed where the combinations of model grid 
spacing (A,„), land-surface properties (L)- and orography (O)- resolution are used 
in a way that addresses their relative importance to convective precipitation. This 
way, the study contributes to the second phase of the Federal Ministry of Educa- 
tion and Research (BMBF) project “HD(CP)* ” which stands for High Definition 
Clouds and Precipitation for advancing Climate Prediction. This project focusses 
on cloud processes by combining current low-resolution global climate simula- 
tions and regional high-resolution simulations. The high-resolution simulations 
facilitate the process-based comparisons with observations, which should give an 
added value to the understanding of cloud processes. The current work is a syn- 
thesis part of the project HD(CP)* . In this part, the land-surface heterogeneity will 
assess and reduce the uncertainties of climate models caused by non-resolved sub- 
surface processes and its feedback to the CBL, the exchange processes between 
CBL and the free troposphere including the cloud development, convection initi- 
ation, and precipitation. In the HD(CP)* framework, the present thesis contributes 
to the investigations about the impact of the land including sub-surface, vegeta- 
tion and anthropogenic structures on the regional climate of central Europe with a 
focus on cloud and precipitation development, intensity, and distribution (Stevens 
et al., 2020). 

This study is not a model evaluation but a model inter-comparison. Nevertheless, 
the assessment concerning the plausibility and the more realistic ICON simula- 
tions are based on qualitative comparison with observational data. 

To summarize the outline: in the next chapter, the literature synthesis of the 
various aspects of model strategies and their limitations are presented (Chapter 2). 
In Chapter 3, the details of the selection of suitable areas and cases are dis- 
cussed. Chapter 4 describes the ICON model system (physics, dynamics and 
their coupling) along with its most interesting and related aspects used in this 
thesis. Chapter 5 presents the details of the ICON model simulations performed 
in this study. It also describes the developed simulation strategy and the traject- 
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ory model LAGRANTO. The results from these simulations are split into two 
chapters: Chapter 6 discusses the first aim of this thesis which is to investigate 
the impact of model grid spacing and land-surface heterogeneity on clouds and 
precipitation. Chapter 7 addresses the identified important triggering mechanisms 
at the most appropriate grid spacing from all the cases and areas. Ultimately, the 
findings are summarised together with an overall conclusion in Chapter 8. 


2. Literature synthesis 


This chapter provides a literature synthesis of the studies analysing clouds and 
precipitation at different grid scales, especially in the numerical simulations. At- 
mospheric processes occur at different length scales. Even with the advent of 
advanced numerical models with much improved representation of physical pro- 
cesses, it is a challenge to decide the optimum model grid spacing, land-surface 
resolution, and suitable parameterization schemes which solve the closure prob- 
lem of the equation of motion entirely efficiently. The atmospheric and land- 
surface based processes play crucial roles in the evolution of convection. There- 
fore, modelling and understanding the processes of clouds and precipitation re- 
quire the selection of a suitable grid scale. According to the scale definition of Or- 
lanski (1975), the relevant scales in this study belong to the meso-ß, meso-y, and 
micro-d, 1.e., the atmospheric processes with a characteristic horizontal scale of ~ 
20 km, 2-20 km, and 200 m - 2 km, respectively whose time scale ranges between 
a diurnal period to few hours (Figure 2.1). 

Further, Inoue et al. (2008) uses the global cloud-resolving model (Nonhydrostatic 
icosahedral atmospheric model (NICAM)) and resolve clouds at a grid spacings 
of 3.5 km and 7 km. A thorough comparison of cloud properties against Japanese 
geostationary meteorological satellite (MTSAT-IR) proves that the overall cloud 
size distribution simulated in 3.5 km grid spacing is closer to the observation in 
contrast to 7 km. However, 3.5 km underestimates the number of small clouds 
and is incapable of producing very large clouds. Moreover, these processes are 
parameterized in the conventional climate models, which is examined as a major 
source of the model errors and uncertainty in the climate projections (Henderson- 
Sellers et al., 1993; Déqué et al., 2007). The approach of explicitly simulating 
the crucial processes at hecto- or kilometre scale is already in practice by the op- 
erational institutes and the scientific community. For example, the studies using 
Convection-Permitting Models (CPMs) like Prein et al. (2015) prove that the ho- 
rizontal grid spacing < 4 km not only outperforms the traditional Large-Scale 
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Figure 2.1.: Scale definition and the characteristic time and horizontal scales of different 
processes (adapted after Orlanski (1975)) 


Models (LSMs) (horizontal grid spacing > 10 km) and present more realistic in- 
sights on climate and regional scales. The modes of explicitly resolving convec- 
tion at a kilometre and finer-scale were already adapted more than twenty years 
back (Weisman et al., 1997). Allin all, the atmospheric mesoscale models are now 
actively used across an extended range of spatial scales with tens of kilometre to 
< 1 km (Chow et al., 2006; Leung et al., 2006; Knote et al., 2010). But even at 
hecto- or kilometre scale the processes like turbulence, cloud microphysics are 
not resolved. Therefore, they have to be parameterized (Bryan et al., 2003). In 
this regards LES runs have been examined for studying clouds (Mellado et al., 
2018). Moreover, the question is how fine grid spacing is enough to represent 
above small-scale processes, how are they represented in the models, and how 
much uncertainty can be expected over a range of grid spacings in LES. 

Several studies are comparing the behaviour of cloud properties across deca- 
, hecto- to a few kilometre scales and provide possible limitations over these 
scales. Chow et al. (2019) addresses the importance of scale-aware parameter- 
ization schemes to study the cross-scale interactions occurring in complex terrain 
within the grey zone continuum (convection, turbulence and topography) at a scale 
ranging from ~ 10 km to ~ 10 m. The grey zone, which 15 also referred as terra 
incognita, corresponds to the grid resolutions where certain features are partly re- 
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solved and subgrid (Wyngaard, 2004). The scale-aware parameterization enables 
a consistent representation of subgrid (turbulence) and grid-scale processes with 
the refinement of grid spacings. For example, a smooth transition from NWP to 
LES-type turbulence closure will be enabled in the so-called scale-aware para- 
meterization schemes (Chow et al., 2019). Exclusively, the uncertainties 10 these 
models may come from initial and lateral boundary conditions, incomplete set of 
parameterized physical processes, or the numerical methods themselves, which 
influences the Convection Initiation (CI), onset and distribution of clouds and pre- 
cipitation (Barthlott et al., 2011). 

Another crucial component is the representation of land surface in these models. 
The main land-surface characteristics can be categorised as land-surface proper- 
ties like vegetation, soil type, soil moisture, and orography. The land surface is 
an integral component to represent in the atmospheric models because of their in- 
tensive impact at different length scales. Land-surface properties and orography 
directly influence the surface temperature and moisture distribution and therefore, 
the partitioning of available energy into sensible and latent heat flux. The het- 
erogeneity of land surface and its associated energy exchange results into a dif- 
ferential heating of the CBL and may modify the atmospheric state by inducing 
secondary circulations which further influences the turbulent transport of heat and 
energy in the CBL and the free troposphere (Avissar and Liu, 1996; Taylor et al., 
2007). That means the clouds and convection are interactively coupled with the 
underlying land surface, which also impacts the initiation and evolution of con- 
vective systems (Banta, 1990; Weckwerth, 2000). The studies based on in-situ 
observations and NWP runs address the importance of land surface-based trig- 
gering mechanisms in the modulation of atmospheric systems (Entekhabi et al., 
1996; Taylor et al., 2007; Garcia-Carreras et al., 2011; Wulfmeyer et al., 2011; 
Khodayar et al., 2013). The land surface-induced thermal and radiant energy ex- 
change processes which occur when there is differential heating caused by land- 
surface heterogeneity, for example, thermally-driven wind system in complex ter- 
rain (Zardi and Whiteman, 2013). Over flat terrain, the impact of variability of 
land-surface properties on convection and turbulence is intensively studied us- 
ing observations (Pielke Sr, 2001; Maurer et al., 2016) and NWP and LES mod- 
els (Shao et al., 2001; Adler et al., 2011; Khodayar and Schädler, 2013; Maurer 
et al., 2016; Liu et al., 2017). The horizontal extent of land-surface anomalies 
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is an important factor for the development of secondary circulation systems. For 
example, the soil-moisture anomalies having a horizontal length scale in the order 
of 2.5 to 10 km tend to generate thermally-driven circulations, for example, the 
atmospheric flows at meso-D and meso-y scales (Figure 2.1) (Shuttleworth, 1991; 
Taylor et al., 2007; Gantner and Kalthoff, 2010). As the interactions between the 
soil and overlying atmosphere are crucial to the climate and weather system over 
a spectrum of spatial and temporal scales (Pielke Sr, 2001; Koster et al., 2004), 
the initial state of soil moisture fields must be estimated correctly (Van Weverberg 
et al., 2010; Barthlott and Kalthoff, 2011). The positive soil moisture-precipitation 
feedback has been consistently proven in observation and model studies (Findell 
and Eltahir, 1997; Schár et al., 1999; Pal and Eltahir, 2001). The wide range of pre- 
cipitation response over the complex terrain has been systematically studied using 
high-resolution simulations (Schneider et al., 2018; Baur, 2019). The model res- 
olution up to 500 m allows to capture thermally-induced circulations occurring at 
micro-Q to meso-y scales. Especially under fair weather conditions, the onset time 
and precipitation amount shows considerable improvement in the high-resolution 
simulations (Schneider et al., 2018). Moreover, the field campaigns like Interna- 
tional H2O Project (IHOP) 2002 (Weckwerth, 2000), Convective Storm Initiation 
Project (CSIP) (Browning et al., 2007), Convective and Orographically-induced 
Precipitation Study (COPS) (Kottmeier et al., 2008; Wulfmeyer et al., 2011) were 
conducted especially to understand the line of processes leading to convective pre- 
cipitation in complex terrain. The low-level convergence zones, heat and moisture 
distribution show significant impact on CBL resulting in deep convection (Wilson 
and Schreiber, 1986; Kalthoff et al., 2009). 

Thus LES is well-suited option to study these small scale processes from micro- 
scale turbulence to the organised mesoscale structures specifically dependent upon 
atmospheric stability, land-surface heterogeneity and orography (Deardorff, 1972; 
Courault et al., 2007; Gantner et al., 2017; Imamovic et al., 2017). On one hand, 
the diversity of the small-scale processes leads to the complexity of their mutual 
interaction, on the other hand, the inclusion of corresponding optimised paramet- 
erization schemes are computationally expensive. For example, in the case of tur- 
bulence parameterization, the ID TKE-based NWP and 3D classical Smagorinsky 
LES comprehend differently as per the applicative area of study. The differences 
on the impacts of convection due to the different turbulence closure are analysed 
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in several studies, for example, Simon et al. (2019). This study shows that the on- 
set time of convection between 1D TKE and 3D Smagorinsky differs almost by 2 
hours. Their research also introduced the Dynamic Reconstruction Model (DRM), 
a mixed model using explicit filtering and reconstruction technique, a turbulence 
closure which results into a shorter extent of the grey zone. 

The intriguing scenario is that over heterogeneous land surfaces and for a given 
model grid spacing, the representation of land-surface heterogeneity is less ap- 
prehended which leads to inaccuracies like the aggregation effect and the dynam- 
ical effect (Giorgi and Avissar, 1997; Mahfouf et al., 1987; Avissar and Schmidt, 
1998). The aggregation effect, which is the result of highly nonlinear subgrid- 
scale processes, primarily affects the surface sensible and latent heat fluxes, the 
dynamics of soil moisture, intercepted water and runoff. The dynamical effect 
results when the heterogeneity induced circulations are not explicitly resolved 
by a model. It influences the CBL structures, evolution of clouds, precipitation 
occurrence and vertical exchange of energy, heat and moisture up to the free tro- 
posphere. Shao et al. (2001) studied the behaviour of these effects by performing 
numerical experiments using the flow model FOOT3D (Briicher, 1997) at atmo- 
spheric resolution of 1 - 4 km in combination with varying soil model resolution 
(1 - 4 km). The results show that at a given atmospheric model resolution, the 
higher land-surface resolution leads to improved or more variable representation 
of surface energy and momentum fluxes and vice versa. That means the land- 
surface induced subgrid variations in surface energy fluxes needs to be taken into 
account in subgrid closure schemes. Furthermore, averaging these local variations 
of the land cover has a much smaller impact on the surface fluxes over homogen- 
eous land surfaces in comparison to that over heterogeneous surfaces due to the 
nonlinear nature of processes and the extended-spectrum of scale variations over 
heterogeneous surfaces (Heinemann and Kerschgens, 2005). 

Moreover, the present study uses ICON-LES, which has the unstructured 
horizontal grid structure and an added value to study the multi-scale interactions. 
This study discusses the land surface-induced circulations and compares the im- 
pact of land-surface heterogeneity and the resolution of the atmosphere on clouds 
and precipitation. In long term perspective, the incorporation of an optimum land- 
surface resolution in climate models will assess and reduce the uncertainty caused 
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2 Literature synthesis 


by subgrid subsurface processes, and its mutual relation with CBL, CI, cloud de- 
velopment, and precipitation. 
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3. Selection of suitable areas and cases for 
convective precipitation 


This chapter describes the approach of selecting suitable areas and the meteoro- 
logical case studies formulated in this thesis. Lightning in thunderstorms can be 
taken as an equivalent to deep convection (Mansell et al., 2007; Leary and Ritchie, 
2009). It gives the general estimate of the occurrence and temporal evolution of 
convective systems. Lightning is defined as the discharge, which could occur as 
cloud-to-cloud (IC), cloud-to-ground (CG), and cloud-to-air (CA). In general, CG 
lightning has a significant impact on our daily lives, in terms of life and property 
loss, power outages and also natural calamities like forest fires. 


Lightning 


The lightning data used here are obtained from the Siemens BLIDS detection 
system, which is a part of the EUCLID network (Schulz et al., 2016). EUCLID 
is a collaboration among the national lightning detecting networks in Europe. It 
aims to identify and detect lightning by using 148 sensors throughout Europe. 
The application of the EUCLID product is not only limited to the meteorology 
but actively in insurance, hydrology, communications, forestry, defence, aviation, 
and hydrology. According to Pohjola and Mäkelä (2013), the implementation 
of a low-frequency (LF) (30 - 300 kHz) lightning detection system results into 
a significantly lower detection efficiency of cloud-to-cloud (IC) lightning. The 
detection efficiency has been confirmed to be 96% for flashes with a peak current 
of 2 kA while the CG lightning strikes can be detected with an accuracy of ~100 
m (Schulz et al., 2016; Piper and Kunz, 2017). Therefore, only cloud-to-ground 
(CG) lightning has been considered in further calculations, and both polarity and 
current have not been taken into account (Kunz et al., 2018). 

To identify the areas with hotspots of Deep Convection (DC), the light- 
ning density during the summer period (MJJAS) of 2000-2018 is investigated. 
Lightning density is defined as the mean daily flash total within a specified grid 
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3 Selection of suitable areas and cases for convective precipitation 


per 25km? 
55°N- | 
1350 
54°N | 
1200 
53°N H- n 
n 900 
US 750 
600 
50°N 
450 
49°N - 
300 
A7°N i 


6°E 8*E 10°F 12°E 14*E 


Figure 3.1.: Total lightning density over Germany during the summer period (MJJAS) of 
2000 - 2018. 


cell (Piper and Kunz, 2017). Figure 3.1 shows Germany wide lightning density 
per 25 km? during the summer period (MJJAS) of 2000 to 2018. 

The higher lightning density over the mountainous or complex terrain like the 
Black Forest, Alps mountains, and Thuringian forest 1s considerably higher in 
comparison to the surrounding flat region. This may be attributed to the orograph- 
ically induced deep convective systems. Additionally, there are smaller areas near 
the southwestern part of Berlin and Harz mountains with the strong number dens- 
ity of lightning strikes. Based on the lightning density (Figure 3.1) and the differ- 
ent types of orographic complexity scale (Figure 6.1), three distinct geographical 
areas have been chosen, namely: Al: the flat terrain near Berlin, A2: the isolated 
mountain range called as Harz mountains with the peak height of 1140 m, and A3: 
the Black Forest mountain range. 

As locally triggered convection is of major interest for this study, two criteria are 
applied to detect the suitable days, that 1s (1) having a weak large-scale synop- 
tic forcing and (i1) a considerable number of lightning strikes over the respective 


1 


areas. As a proxy for weak synoptic forcing, we used a threshold of «10 ms” at 
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3 Selection of suitable areas and cases for convective precipitation 


850 hPa. For the selection of suitable days, the number of lightning strikes and the 
horizontal wind at 850 hPa from ECMWE Re-Analysis (ERA)-Interim (Berrisford 
et al., 2011; Dee et al., 2011) dataset are used. Figure 3.2 shows a typical example 
of the outcome of the above criteria calculated over Al where the set of markers 
represent the different years, and the numbers are the days of the summer period 
(MJJAS) showing considerable lightning density. Similarly, using this criterion, 


Horizontal wind at 850 hPa 


100 200 300 400 500 600 
Number of lightning strikes 


Figure 3.2.: Distribution of summer days (MJJAS) over Al in the period 2012 - 2017 on 
the basis of areal mean horizontal wind speed at (at 850 hPa) and the number 
of lightning strikes. The number indicates the days in the summer period (e.g. 
1 indicates 01 May). Number 87 represents one of the selected case, dated 26 
July 2012. 


two suitable cases are chosen for each area, that means, in a total of six cases. To 
get an impression of the observed precipitation patterns on the investigated days, 
the Radar Online Adjustment (RADOLAN) RW composites are shown in the fol- 
lowing section. However, notably, this study focuses on the inter-comparison of 
simulated precipitation for different model grid spacings and land-surface res- 
olution, not on their validation against the observation. As the name suggests, 
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3 Selection of suitable areas and cases for convective precipitation 


RADOLAN from the DWD is a forum which provides Radar-based quantitat- 
ive precipitation estimation products which are of exceptionally high temporal (5 
minutes) and spatial resolution (1 x 1 km?) covering all hydrological catchment 
areas of Germany. The high spatio-temporal resolution provides data input to 
flood risk management and is used for other hydro-meteorological and climatic 
applications (Bartels et al., 2004). The precipitation products in RADOLAN RW 
products have been used in this study. These products are based on the combin- 
ation of quantitative measurements with C-band Doppler radar and hourly rain 
gauge measurements. Here, RW denotes hourly adjusted radar data. 

In the following section, the general synoptic situation of one example case 
from each area is briefly discussed using the composites of the observed daily- 
lightning density and accumulated precipitation. The other cases are shown in the 
appendix (Figure B.2, B.3, B.4). 


3.1. A1: Flat terrain 


Case: 26 July 2012 


The spatial distribution of daily lightning density per 25 km? and accumulated 
precipitation (in mm) from RADOLAN RW on 26 July 2012 over Al are shown 
in Figures 3.3a and 3.3b. This day has scattered but intensive precipitating cells 
restricted mainly in the east and southeast of Berlin. The high lightning density 
indicates that the precipitation result of deep convection. 


3.2. A2: Isolated orography 


Case: 09 June 2018 


This day is characterized by a locally induced precipitating system over A2 (Harz 
mountains) (Figure 3.4a and 3.4b). There are mainly two sets of convective sys- 
tems observed over this region. The first system was formed along the Harz 
mountains around 1100 UTC and the second subsequent system was formed in 
the southwestern part of Harz mountains around 1300 UTC. Both of the systems 
show intensive lightning density indicating deep convective activity. 
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3.3. A3: Complex terrain 
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Figure 3.3.: Horizontal distribution of lightning density (a) and observed (RADOLAN 
RW) daily accumulated precipitation (b) over A1 on 26 July 2012. 
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Figure 3.4.: Same as Figure 3.3 but over A2 on 09 June 2018. 
3.3. A3: Complex terrain 
Case: 12 Aug 2015 


The composites of daily lightning density per 25 km* and accumulated precip- 
itation (in mm) from RADOLAN RW on 12 Aug 2015 over A3 are shown in 
(Figure 3.5a and 3.5b). It shows the intensive precipitating cells formed over the 
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3 Selection of suitable areas and cases for convective precipitation 
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Figure 3.5.: Same as Figure 3.3 but over A3 on 12 Aug 2015. 


peaks of the Black Forest, especially over the southern Black Forest, the precipit- 
ating cells are convective and intense. 
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4. ICOsahedral Non-hydrostatic (ICON) Model: 
description of its key features 


This chapter describes the general features of the ICON model. The ICON is 
a nonhydrostatic fully compressible general circulation model that has been de- 
veloped at German Weather Services (DWD) and MPI-M (Zängl et al., 2015). 
ICON is used at DWD operationally as a global numerical weather prediction 
model since January 2015. Different communities for the assessment of the un- 
certainty in the representation of mesoscale circulations in comparison to other 
NWP models (Heinze et al., 2017). The authors also aim at high resolution, re- 
gional LES and idealised (Dipankar et al., 2015; Silvers et al., 2016) simulations 
to study the subgrid forcing of the physical processes, for example, turbulence, 
cloud-based processes. ICON also intends to a better representation of the ex- 
change processes between different air masses along with an improved represent- 
ation of the land surface and subgrid-scale heterogeneities. Depending upon the 
applications and scale study, ICON is further designed in three basic packages; 
namely, 1) Climate predictions (ICON-Global), 11) Numerical weather prediction 
(ICON-NWP) and iii) Large-eddy-simulations, (ICON-LES), which is also called 
as ICOsahedral Nonhydrostatic - Large Eddy Model (ICON-LEM). 

ICON-LEM uses physics designed for LES ın Limited Area Mode (LAM) 
usually at the scale of O(100m). The LES physics package was introduced in 
ICON (Dipankar et al., 2015) within the framework of HD(CP)* which stands for 
High Definition Clouds and Precipitation for Climate Prediction project. HD(CP) 
aimed at the improvement of representation of clouds, convection and precipita- 
tion processes in the climatic prediction by utilising very high O(100m) resolution 
simulations. These simulations provide a better understanding of embedded and 
complicated resolved processes which are parametrized in the coarser model sim- 
ulations. Furthermore, ICON-LEM model simulations also provide lateral bound- 
ary conditions for other limited area forecasts which are used in the nested simu- 
lations for a better representation of boundary and environmental flows. 
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4 ICOsahedral Non-hydrostatic (ICON) Model: description of its key features 


ICON follows the sets of equations of prognostic variables suggested by Gass- 
mann and Herzog (2008) and the two dimensional Lamb transformations which 
transform the nonlinear horizontal momentum advection into a vector invariant 
form as explained in Zängl et al. (2015). These set of equations describe a two- 
component system which consists of dry air and all the three phases of water. The 
prognostic equations are solved for horizontal velocity components normal to the 
triangular edges v„, the vertical wind component w, the total density of the air 
mixture p, virtual potential temperature 0, and the specific masses and number 
densities of tracers q; (i = 1,2,3,.., N;, where N; is the total number of traces). 
According to Zängl et al. (2015), the basic equations system used in ICON are as 


follows: 

3: t3. 6 + fy twa = cpa ys Fn); (4.1) 

Ow Ow oll 
E హై M = N), 8 (4.2) 

ð 
F +V (vp) = 0, (4.3) 
d 0, R 

E + V.(vp8,) = Ó. (4.4) 


For a simplified numerical treatment of the terms representing vertical sound wave 
propagation, the equation 4.4 can be also written as: 


on R, I 


Ar erga V 009) =Ó (4.5) 


where Kp = (vn? + vj?) is the horizontal component of kinetic energy per mass 
unit, v; is the reconstructed tangential velocity component, ¢ is the vertical vorti- 
city component, f is Coriolis parameter, g is the acceleration due to gravity, F (vn) 
is the source term for horizontal momentum, v, is horizontal velocity compon- 
ent, v is the full three dimensional wind vector, II is the Exner function, Ry is 
the gas constant of dry air, Cpa and 64 are specific heat capacities of dry air at 
constant pressure and volume, respectively, and Ó is the diabatic heat source term 
(Wan et al., 2013; Zängl et al., 2015). Moreover, the derivative, 2 denotes a ho- 
rizontal derivative in edge-normal direction. In ICON, local mass conservation is 
achieved through flux form of equation 4.3 and 4.4. For time integration, a two- 
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4.1. Horizontal grid 


time level predictor-corrector scheme is used, which is explicitly done except for 
the terms describing sound-wave propagation (for details see Zängl et al. (2015)). 
To achieve a mass-consistent tracer transport, the air mass fluxes are aggregated 
over the small-time steps in the dynamical core and then passed to the transport 
scheme. This tracer transport 1s done using a flux-form semi-Lagrangian scheme 
(Miura, 2007). 

The ICON grid system comprises of horizontal grid (as described in Sec- 
tion 4.1) and the vertical grid setup (see Section 4.2). According to Linardakis 
et al. (2011) the horizontal grid discretises the sphere surface on a triangular or 
hexagonal grid, and the vertical grid is discretised using a set of horizontal layers 
along the sphere radius. 


4.1. Horizontal grid 


The horizontal grid of ICON is one of the most featured advantages of this model. 
The first step to create the horizontal grid is projecting a regular icosahedron onto 
the sphere in such a way that two vertices fall onto North and South Poles (Figure 
4.1). Figure 4.1 shows an example of the regional limited ICON grid where the 
light blue grid indicates the icosahedron defined base grid, which comprises of 
20 triangular faces, 13 edges and 12 vertices. The first refinement of this root 
grid is the dark blue grid where the edge centres form the new vertices. There is 
another refinement done only for the northern hemisphere (green grid) and based 
on this grid, another mesh refinement is done over the European region (red grid). 

The desirable resolution is achieved by triangulating the triangular faces of the 
icosahedron (cf. Figure 4.2a), which is called Delaunay-Voronoi triangulation 
method, which is nothing but a repeated subdivision of the triangular cells of a 
spherical icosahedron (cf. Figure 4.2b) into smaller cells. This results into the 
division of each great circle arc of a projected icosahedron face into n arcs of 
equal length (cf. Figure 4.2c) and each icosahedron face into n,^ small triangles. 
For details, see Sadourny et al. (1968). This process is called root division, and 
the resulting mesh is grid level 0. With the further mesh refinement (by bisecting 
each spherical triangle edge, k) and connecting the midpoints by great circle arcs, 
it leads to ‘RnBk’ grid (Figure 4.2). Figure 4.2d,e show the two examples of R2B0 
and R2B2 grids. 
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4 ICOsahedral Non-hydrostatic (ICON) Model: description of its key features 


Figure 4.1.: The icosahedral grid of the ICON model. The grid level 0 or icosahedron base 
grid (light blue lines) has 20 triangular faces, 30 edges and 12 vertices. The 
first horizontal grid refinement is shown by dark blue lines. The green grid 
shows the refinement done over the northern hemisphere, and the third one is 
done over Europe, shown by red lines (Source: Giorgetta et al. (2017)). 


(b) (c) (d) 


(f) 


Figure 4.2.: Construction of ICON horizontal grid. The icosahedron (a) is constructed on 
a sphere (b). The triangle edges are bisected into n equal sections. The new 
edge points are connected great circle arcs to generate n? spherical triangles 
within the root triangle (c). After further mesh refinement, the aimed grid is 
obtained (e). (f) compares the occurrence of polar singularities of latitude- 
longitude grids against the complete avoidance of it in (e) for the ICON grid 
construction (Source Linardakis et al. (2011)). 


Figure 4.2e highlights the advantage of triangular grid in ICON. It removes the 
polar singularities of latitude-longitude grids and allows a consistent uniformity 
in resolution over the globe (for further advantages of the triangular grid, see 
Gassmann (2011)). 

There is an interesting distinction and definition of the term ‘resolution’ in ICON. 
According to Heinze et al. (2017), the resolution in ICON is defined as the square 
root of the average cell area of the icosahedral grid of ICON, which is approxim- 
ately 1.5 times the resolution on a corresponding regular grid. This means that for 
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4.2. Vertical grid 


the resolution R,B;, the total number of cells (ne) , edges (ne) and vertices (n,) 
will be: 


Ne = 20n^4*, n, = 30n74*, n, = 10n^4* +2. (4.6) 
Further, the effective mesh size Ax is: 


Ax = Va; = fan, (4.7) 
where a. is the average cell area and r, is the radius of the Earth (Zängl et al., 
2015). Altogether, the definition of effective mesh size or grid resolution could 
be also based on the square root of the area of the triangle, the length between 
the two cell-centres, or the edge length of the cell. According to Dipankar et al. 
(2015), the approximate relation between the triangle edge length (A;) and the grid 
resolution (A) following the above definitions will be: A; = 0.67A;, A» = 0.58A;, 
A3 = A, respectively, where A4 > A4 > Ad. 

But there are also slightly different definitions of "resolution" in ICON used by 
other researchers. According to Wan et al. (2013), in the grid level 0, there are only 
12 icosahedrons which are surrounded by five triangles, also called as pentagon 
points, the rests are surrounded by six triangular cells. This grid distribution res- 
ults in inequality in the cell areas and the edges of triangles. Quantitatively the 


reason is that the average cell centre distances are smaller by a factor of about 
4 
3 
bours. For the present study, the grid resolution is taken as the triangle edge length 


on a triangular grid because every grid cell point has only three nearest neigh- 


which is termed as ' model grid spacing’ (Aj) (this should not be confused with 
the model resolution). Furthermore, C-staggering is applied to the triangular cells 
by putting mass and temperature at the circumcenters of triangles (Bonaventura 
and Ringler, 2005; Wan et al., 2013; Giorgetta et al., 2015; Linardakis et al., 2011) 
(Figure 4.3). 


4.2. Vertical grid 


The vertical grid of ICON comprises of a set of layers which possess the hori- 
zontal two-dimensional grid structure, shown in Figure 4.4. It employs a Lorenz- 
type staggering where the vertical velocity is defined at the boundaries of layers 
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Figure 4.3.: Positioning of prognostic variables on ICON horizontal grid with primal cell 


(triangular) and dual cell (hexagonal). The dual edges are bisect and ortho- 
gonal to the primal edges (Source: Linardakis et al. (2011)). 


(half levels). At the same time, other prognostic variables are placed in the centre 
of the layers (full levels). ICON has a height-based terrain-following hybrid struc- 
tured vertical coordinate system and the smooth level vertical coordinate imple- 
mentation (SLEVE) (Leuenberger et al., 2010). The conventional Gal-Chen and 
Somerville (1975) formulation is also available alternatively. The main advant- 
age of this vertical coordinate system is that it allows a faster transition to smooth 
levels in between the upper troposphere and lower stratosphere in comparison to 
the traditional height-based Gal-Chen coordinate system. 

Depending upon the specified model top height, the number of height levels and 
the stretching factor, ICON defines the set of vertical levels. These definitions 
of vertical grid staggering are controlled by the use of the respective namelist 
parameters (number of levels, stretching factor, model top height, minimum layer 
thickness, flat height (the height above which the vertical coordinate surfaces are 
flat) (Namelist Overview, 2015). The stretching factor controls the thickness and 
therefore, the distribution of vertical levels. One example of the vertical level 
distribution of ICON-LEM simulations performed in this study is shown in Figure 
4.5. The detailed model configuration is discussed in Section 5. 


4.3. ICON physics, dynamics and their coupling 


In ICON, the physical processes are categorised as slow- and fast- physics forcing 
depending upon how frequently they are called in the simulation. In other words, 
the processes whose time scale is comparable or relatively shorter (longer) than 
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4.3. ICON physics, dynamics and their coupling 


Figure 4.4.: Illustration of the vertical structure of the primal cells. As shown, the half 
levels (k + 1/2) correspond to n vertical levels (Source: Linardakis et al. 
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Figure 4.5.: Distribution of vertical half levels (blue) and the layer thickness (red) in one 
of the ICON simulations performed over flat terrain (A1). 


the model time step are called as fast (slow)-physics processes. The examples 
of fast physics are saturation adjustment, surface transfer, land/lake/sea-ice based 
parametrization, turbulent diffusion and microphysics. Slow physics are convec- 
tion, cloud cover, radiation, non-orographic gravity wave drag and sub-grid-scale 
orographic drag. Since the slow physics tendencies are called less frequently, they 
are stored to be integrated with the governing equation. In contrast, the fast phys- 
ics are called every model time step so that it updates the prognostic variables 
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4 ICOsahedral Non-hydrostatic (ICON) Model: description of its key features 


sequentially and therefore, they do not need to provide the tendencies terms to 
the model governing equations (Zängl et al., 2015; Dipankar et al., 2015; Heinze 
et al., 2017). Since ICON is designed to solve prognostic equations based on the 
Exner function (II) rather than temperature (T), the temperature tendencies have 
to be converted into tendencies of II first. Also, the moisture tendencies are not 
treated as a forcing term during the tracer advection but in a split manner and then 
further added to the newly updated moisture variables. 


4.4. Land-Soil Model TERRA_ML 


The atmosphere and the underlying land-surface 15 strongly coupled, which is 
represented in numerical models as surface fluxes. In ICON, the soil-vegetation- 
atmosphere-transfer component TERRA_ML (Schrodin and Heise, 2002; Heise 
et al., 2006) is used for the exchange of momentum-, heat-, moisture-, mass- 
fluxes between land-surface and atmosphere. It forms the boundary condition 
for the atmospheric part of ICON. TERRA_ML is the multi-layer version of the 
soil model TERRA. It considers the physical processes like radiation, biophys- 
ical control of evapotranspiration, heat and soil-water transport, snow formation 
and the coupling with the atmosphere based on a multi-layer concept for the soil. 
Usually, most of the characteristic parameters of the soil model, for example, heat 
and water storage capacity, porosity, are strongly coupled with the soil texture. 
With this categorisation in TERRA_ML, it gives eight different soil types. The 
multi-layer version not only makes the layer thickness independent of the respect- 
ive soil type but also avoids the different soil layer designation into the thermal 
and hydrological section of the model (for details see Doms et al. (2011)) 


4.5. ICON-LEM standard setup 


As stated before, ICON-LEM is the LES setup of ICON designed for very high- 
resolution simulation, which is especially done in LAM. Besides having high po- 
tential in a line of improving the understanding of moist convective processes and 
their parameterization in climate models (Heinze et al., 2017), ICON has proven 
its eddy-resolving abilities (Dipankar et al., 2015). However, the implementation 
of a three-dimensional turbulence scheme on a triangular grid on a conservative 
basis is described as a non-trivial step because many interpolation operations have 
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4.5. ICON-LEM standard setup 


to be performed. Additionally, there are new approaches mandatory for the rep- 
resentation of the subgrid-scale turbulence and other complicated microphysical, 
cloud-based processes (Dipankar et al., 2015). The adaptation of the simulations 
using the standard ICON-LEM setup is described in detail in the next chapter 5. 
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5. Adaptation of ICON for sensitivity simulations 


5.1. Domain configuration 


This chapter describes the ICON simulations performed in this study. The ICON 
simulations are performed in LEM mode, with the model grid spacings (4500077. 
Aos500m» A1250m; A625m; A312m and A156m) for a comparison of the impact of differ- 
ent model grid spacing (A;) and land-surface resolution on moist convection. The 
six model grid spacings (ర) are used in a nested setup, starting from the parent 
domain of grid size 5000 m, going to a grid size of 156 m, in the ratio of two. A 
typical example of the domain configuration for a case simulated over flat terrain 
(Al) on 26 July 2012 is shown in Figure 5.1. The setup uses 90 vertical levels 
with the model top height at 20,000 m, and the lowermost minimum layer thick- 
ness as 10 m which results into the stretching factor of model levels that equals to 
0.9 (Figure 4.5). 


5.2. Model physics setup 


There are limited studies which address the suitable resolution when there is 
a transition from 1D turbulence scheme of NWP to 3D turbulence scheme of 
LES (Honnert and Masson, 2014; Cuxart, 2015). According to Honnert and Mas- 
son (2014) the criterion of an optimal horizontal grid spacing (A,) scaling should 
be based on the size of thermal plumes which is equal to the sum of boundary layer 
height h and the depth of cloud layer h,. This study proposed two separate critical 
resolution for free and forced CBL which is Ax € 0.5(h+h.) and Ax € 3(h+ he), 
respectively. However, ICON-LES turbulence scheme (3D Classical Smagorinsky 
scheme) assumes isotropic turbulence which is not a good approximation espe- 
cially for shear-driven turbulence and therefore a scale adaptive turbulent scheme 
is recommended. The choice of physics packages in this study is based on a few 
tests and for a broader range of comparison of the grid spacings (Aj). Therefore, 
the transition in turbulence closure is done for A1250m and the finer grid spacings, 
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5 Adaptation of ICON for sensitivity simulations 


which are referred as LES model grid spacings. This scale covers the atmospheric 
flows ranging from mesoscale y to microscale &. Mesoscale simulations still do 
not explicitly resolve the atmospheric motions (and turbulence) occurring over an 
hourly or lower time range. But the chosen LES model simulation range explicitly 
resolve most of the deep convection and turbulent eddies. The selection of LES 
turbulence for model grid spacing from Aj259 to Aj56m allows us to compare the 
sensitivity of model grid spacing on deep convection over a broader range with a 
simple and uniform setup. 


Asooom 
A2500m 
Ai250m 
Ae25m 
A312m 


A156m 


10°E 11°E 12°E 13°E 14°E 15°E 16°E 17°E 


Figure 5.1.: Model domain for the case over flat terrain, area near Berlin, dated 26 July 
2012. This is a setup of one way nesting starting from Asooom nested down to 


A156m. 


Since these packages are designed for studying the eddy-resolving as well as also 
subgrid-scale processes, there are two sets of processes which have to be expli- 
citly and implicitly represented in the model simulation. Firstly, the explicit phys- 
ical processes are categorised into the following physics packages. For example, 
convective schemes, subgrid-scale orographic effects (this controls blocking and 
gravity wave drag formation), non-orographic gravity wave drag. The second cat- 
egory is for the processes, which need to be parametrized, and which are called 
under the following physics packages: 


- Cloud microphysics scheme: Double-moment microphysics scheme based 
on Seifert and Beheng (2001) 


- Subgrid-scale turbulence scheme: 
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5.3. External parameter: EXTPAR 


e Prognostic Turbulent Kinetic Energy (TKE) COnsortium for small-scale 
MOdelling (COSMO) (for NWP) (Raschendorfer, 2001) 


e 3D Smagorinsky diffusion scheme (for LES) based on Smagorinsky 
(1963) with the modifications by Lilly (Lilly, 1962) in order to consider 
the thermal stratification (buoyancy along-with wind shear) 


- Cloud cover scheme: 


e Diagnostic Probability Density Function (PDF) cloud cover scheme (for 
NWP) 


e ‘All-or-nothing’ which is also called a grid-scale cloud cover scheme 
(for LES) 


- Land-soil model: TERRA-ML (Schrodin and Heise, 2002). It consists of 
eight irregularly spaced soil levels and eight soil types. In the current model 
setup, ntiles = | has been used, which means that the tile approach is turned 
off and only the dominating land-surface type of a grid cell has been taken 
into account. 


- The ICON simulations performed in this study do not use a lake or sea-ice 
model. 


The Rapid Radiation Transfer Model (RRTM) scheme is applied for radiation 
(Mlawer et al., 1995). TERRA_ML is used as multi-layer surface scheme (Doms 
et al., 2011). This is the version without subgrid land-cover variability, and in 
this configuration, the vertical interpolation of soil variables are needed. Land- 
use data, orography and soil type are taken from the high-resolution observational 
dataset, which is also interpolated on respective ICON grid. The details of these 
external parameters are given in Section 5.3. 


5.3. External parameter: EXTPAR 


At the step of initialising the surface fields, NWP models need to generate the 
geographically localised dataset for example terrain height, plant cover, land use. 
External Parameter for Numerical Weather Prediction and Climate Application 
(EXTPAR) is a software tool used in COSMO and ICON to perform this task. 
It takes the raw input dataset and interpolates it on the target grid (here, ICON). 
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Parameter Dataset Source Resolution 
Globcover 2009 European Space T 
ns (GlobCover, 2009) Agency (ESA) en) 


Ministry of Economy, 
Trade, and Industry of 
Japan (METI) / National 


Advanced Spaceborne 
Thermal Emission and 


Ones apy Ra aes Aeronautics and Space Pa) 
(ASTER, 2011) Administration, United 
States (NASA) 
Harmonized World Soil 
Soil-type Database (HWSD) a 30” (900 m) 


(Fischer et al., 2008) 


Table 5.1.: Input raw datasets for ICON 


In this way, the different raw datasets are aggregated to the target grid, by con- 
sidering all the raw data elements, which lie within the target grid element. For 
details see Asensio et al. (2018). The source-raw datasets of some of the external 
parameters and their resolution is tabulated in Table 5.1. 


5.4. Initial and boundary conditions 


The initial and lateral boundary conditions for ICON simulations are taken from 
different types of models. Out of the six cases, one case over flat terrain (A1) 
26 July 2012 is initialised using the surface fields from COSMO-DE and at- 
mospheric fields from ECMWF-Integrated Forecasting System (ECMWF-IFS). 
COSMO-DE is a NWP model having a spatial resolution of 2.8 km and 50 vertical 
levels. Here, DE refers to its domain setup over Germany, Switzerland, Austria 
and parts of some other neighbouring countries (for details see Deutscher Wet- 
terdienst (DWD) (2016)). The ECMWF-IFS model operational analysis dataset 
is with a horizontal resolution of ~9 km. With this, the ICON-LEM is relaxed 
towards six-hourly ECMWF-IFS analysis in a 20 km wide nudging zone. The 
nudging zone width of the nested domains is 8 grid points wide (cf. Heinze et al. 
(2017)). 

The other cases use the operational ICON analysis product ICON-Europe 
(ICON-EU) for initial and boundary conditions having a resolution of 6.5 km and 
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60 vertical layers. The boundary of the parent domain (Asooom) 15 relaxed towards 
three hourly ICON-EU assimilation forecast products. ICON-EU products are 
available to registered users via DWD’s PArallel MOdel data REtrieve from the 
Oracle Databases (PAMORE) service (for details see DWD-PAMORE (2015)) 

The above-listed model products need to be interpolated on the correspond- 
ing ICON grid in order to be used as initial and boundary conditions. For this 
purpose ICON-tool from DWD is used (for details see Prill (2014)). For example 
ICON-EU data has been interpolated on the respective ICON-LEM grids by using 
a Radial Basis Function (RBF) interpolation method (Peixoto and Barros, 2014). 

These six different days are simulated each for 24 h periods (00-24 Co- 
ordinated Universal Time (UTC)). A summary of the days, with some important 
notations and differences, is shown in Table 5.2. 


Al A2 A3 
Areas Å ; 
Flat terrain Isolated mountain range | Complex orography 
26.07.2012 09.06.2018 12.08.2015 
oo 17.09.2017 17.09.2017 29.05.2017 
Model resolution 5000 m, 2500 m, 1250 m, 625 m, 312 m, 156 m 
Duration 00-24 UTC 00-24 UTC 00-24 UTC 
COSMO-DE and ECMWF-IFS ICON-EU ICON-EU 
Inca poudion ICON-EU ICON-EU ICON-EU 
ECMWF-IFS ICON-EU ICON-EU 
en ICON-EU ICON-EU ICON-EU 


Table 5.2.: Outline of simulated cases 


5.5. Simulation strategy 


The general structure of the set of simulations is discussed in section 5.1, 5.3 
and 5.4. This part discusses the simulation strategy of all the cases, along with 
their sensitivity simulations (see Figure 5.2). In all the cases, A156m of the set 
of control runs is taken as the reference run. The control run simulations have 
been performed for all six cases. Along the control run line, the model grid spa- 
cing (blue) and land-surface resolution (green) have been changed simultaneously 
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(shown along the diagonal line, Figure 5.2). Furthermore, the sensitivity simula- 
tions have been performed by keeping the land-surface resolution (<res>: 5000 
m and 1250 m) constant and only changing the model grid spacing. Here, land- 
surface parameters can be broadly categorised into groups. The first group com- 
prises of parameters related to land-surface properties (denoted as L) for example 
plant cover, land use, the fraction of land and lake, roughness length. The second 
group consists of the parameters related to orography (denoted as O) for example 
terrain height, slope angle, slope azimuth. This way, the sensitivity simulations 
have been done further in three ways, firstly the land surface properties and oro- 
graphy parameters are modified altogether (denoted as LO_<res>). Secondly, the 
land-surface properties are only modified, but the resolution of orography-related 
parameters are the same as the corresponding model grid spacing (denoted as 
L_<res>). Thirdly, the orography-related parameters are modified, but the resolu- 
tion of land-surface properties are used as the corresponding model grid spacing 
(denoted as O_<res>). 


Model grid spacing (m) 


Land-surface resolution (m) 


Figure 5.2.: ICON-LEM simulation strategy. Terms in blue colour represent model grid 
spacing (Aj) and in green colour represent land-surface resolutions (LSR). 
Ref (in red) which is 156 m control run, denotes the reference run. 


That means the land-surface properties and orographic parameters are coarsened 
sequentially while the atmospheric variables are used at high resolution. This 
modification is done by interpolating the land surface parameters from 5000 m to 
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the finer resolution 2500 m, 1250 m, 625 m, 312 m, 156 m in the case of LO_5000 
m, L_5000 m and O_5000 m and from 1250 m to 625 m, 312 m, 156 m in the 
case of LO_1250 m, L_1250 m and O 1250 m. This modification is done by using 
the DWD-ICON-tools. The interpolation is done using the option of Barycentric 
interpolation algorithm. The main advantage of this interpolation method is that 
it interpolates the data from triangular grid exactly on the target triangular grid. 
Therefore there is no over- and undershoots of the resulting interpolated paramet- 
ers, which happened to be the case with other interpolation methods (RBF and 
Nearest-Neighbour method) (for details see Prill (2014)). 


5.6. LAGRANTO: LAGRangian ANalysis TOol 


LAGRANTO stands for LAGRangian ANalysis TOol (Wernli and Davies, 1997). 
It is a tool designed essentially to calculate forward and backward trajectories of 
air parcels to identify the flow structures of air masses. This tool is frequently 
used in the study of extratropical storm tracks, tropical cyclones, warm conveyor 
belt, identification of atmospheric blocking, dispersion of pollutants and in the 
characterisation of other transport regimes (Sprenger and Wernli, 2015; Steinfeld 
and Pfahl, 2019; Babié et al., 2019). The essential theory of the trajectory calcu- 
lation in LAGRANTO is the same as used in other trajectory tools, for example 
FLEXTRA trajectory model (Stohl, 1999), Hybrid Single Particle Lagrangian In- 
tegrated Trajectory model (HYSPLIT) (Draxler and Hess, 1998) etc. The basic 
trajectory equation, which they solve numerically is: 


Dx 

[m u(x), (5.1) 
where x = (A, Ø, p), u = (u, v, అ. x is the position vector in the respective 
geographic coordinate and u represents 3-D wind field. Though LAGRANTO 
differs in terms of spatial discretisation, treatment of lower boundary and also 
its flexible objective criteria of selection of trajectories. It also allows the user 
to calculate some diagnostic parameters like potential vorticity, diabatic heating 
rate, divergence and potential temperature etc. LAGRANTO is flexible in the 
sense that it allows the user to calculate forward and backward trajectories for 
different global and NWP models (European Centre for Medium-Range Weather 


35 


5 Adaptation of ICON for sensitivity simulations 


Forecasts (ECMWF), COSMO, Weather Research and Forecasting (WRF) and 
UK Met Office Unified Model (MetUM)). 

In this study, LAGRANTO is used to identify the source region of the convective 
system. To achieve that the backward trajectories starting from the cloud base 
or inside the cloud itself are calculated. The adaptation of the LAGRANTO 2.0 
version has been done for the ICON model with the assistance of M. Sprenger 
from ETH Zürich. The version LAGRANTO-ICON works on horizontal regu- 
lar latitude-longitude. Therefore, the required trajectory fields are interpolated 
from ICON to regular latitude-longitude grid. LAGRANTO uses 3-D wind fields 
namely (zonal wind u (in m s~!), meridional wind v (in m s^!) and vertical wind 
Q (in Pa s-)). Additionally, it needs surface pressure (PS) (in hPa) for the con- 
figuration of a level type, which further calculates the full 3-D pressure (P) on 
the model fields and identifies the trajectory points intersecting the underlying 
orography. To study the local triggering, the vertical coordinate in meters is in- 
tended. Therefore, LAGRANTO needs the information of orography and the 3-D 
geopotential height (z mnc) for all grid points. The new version is adapted to work 
also with a very high temporal resolution of model fields (which in this study is 
up to 1 to 10 minutes). With the above adaptations and the required trajectory 
fields, LAGRANTO - ICON 15 used extensively in this study. This flexibility of 
LAGRANTO makes it a novel choice for research doing the calculation of air 
parcel trajectories. For details of the different functionalities and their usage in 
LAGRANTO see Sprenger and Wernli (2015). 
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6. Impact of model grid spacing and land-surface 
resolution on clouds and precipitation 


The first aim of this study is to investigate the impact of model grid spacing (A,,) 
and land-surface resolution on convective precipitation using ICON. In the first 
step, the selected areas are classified based on land-surface heterogeneity using 
the parameters representing orography, soil type, and vegetation (Section 6.1). 
It gives an overview of the impact of model grid spacing (Aj) and land-surface 
resolution on convective precipitation (Section 6.2) analysed for the six simulated 
cases (selected as described in Section 3). Moreover, it also describes the respons- 
ible factors and processes which may cause the differences across LES model grid 


spacings (612507 to 615677). 


6.1. Classification of land-surface heterogeneity 


The criteria for selecting the suitable geographical areas and an overview of the 
cases and their synoptic is discussed in Chapter 3. Figure 6.1 shows the hori- 
zontal snapshot of the orographic structures of these areas located in three differ- 
ent regions of Germany. Considering the different physical characteristics of the 
specified region, the land-surface heterogeneity can be categorised broadly into 
two groups, namely, the land-surface properties, for example, vegetation (Leaf 
Area Index (LAI), plant cover), and the orography (Meißner et al., 2007). Fig- 
ure B.1 shows the spatial pattern of modelled orography, Transpiration Area In- 
dex (TAI) and soil type across the six chosen model grid spacings (4500077 A2500m» 
A1250m; A625m> A312m and 15670) of each area (Al, A2 and A3). The orography 
shown here is a model output (Figure B.1a-c). Notably, the ICON model output 
orography is usually produced by digital filtering with a V*-diffusion operator. 
But in the simulations performed in this study, there is no smoothing applied to 
the orography. That means the original orography interpolated on the respective 
model grid spacing (Aj) is used (Figure B.la-c). Furthermore, TAI (Figure B.1d- 
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DORA RR d, ER Er PPIX PAD DOCE i rd 
6°E 8°E 10°F 12°E 14°E 


Figure 6.1.: Orography map (colour coded) with the selected areas: Al (flat terrain), A2 
(isolated mountain range), and A3 (complex terrain). The thin solid line 
marks the national boundary of Germany. 


f) which is the product of plant cover and maximum LAI is a good represent- 
ative of vegetation. Lastly, the third characteristic parameter shown here is the 
soil type (Figure B.1g-1). There are eight different soil types used in TERRA 
- Multi Layer (TERRA-ML) (Schrodin and Heise, 2002), which are: ice, rock, 
sand, sandy loam, loam, loamy clay, clay, peat (seawater and sea ice) which are 
denoted over a scale from | to 8, respectively (Section 4.4). 

The flat terrain area, Al, in Eastern Germany is mainly characterised by a variety 
of land-surface properties, for example, a set of soil type composition (10.5 % 
sand, 40 % sandy loam, 41.5 % loam, 0.4 % clay, 6 % peat, and 1.7 % water) 
and TAI. The second area, A2, up-rises on the orographic-complexity level and 
possesses a higher elevated mountain range, called Harz mountains and also the 
heterogeneous soil type (52.6 % sand, 0.9 % sandy loam, 45 % loam, 0.8 % loamy 
clay, 0.7 % peat, and 0.3 % water) located close to its highest mountain peak (1132 
m), called Brocken. The third area, A3, which is a complex mountain range, called 
the Black Forest, located in south-west Germany and has the highest peak of 1484 
m (Feldberg). A3 consists of a range of orographic as well as vegetation ranging 
from Rhine valley, parts of Swabian Jura and Black forest mountain ranges. A3 
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also has heterogeneous soil type composition with 18 % sand, 3 % sandy loam, 
71 % loam, 1 96 loamy clay, 6 % clay and 0.5 % water. 

Now, in order to quantify this spatial variability of land-surface heterogen- 
eity across the model grid spacings, the areal mean Coefficient of Variation (CV) 
of three physical parameters (orography, soil type and TAI) (Figure 6.2) of the 
three selected geographical areas is calculated. CV is a measure of relative vari- 
ability of a parameter and can be expressed as the ratio of its standard deviation 
(0) to the mean value (u). Therefore, CV for any parameter (P) over an area 
(4;=1,2,3) at a grid spacing (A,) will be: 


౮౧||. 
PA: „100%, (6.1) 
A 


CV pa; = 


where Ltp|4 is the mean value of a given parameter (P) over all the three areas (A). 
CV is expressed as a percentage value. The higher the CV of the given parameter 
(orography, TAI and soil type), the larger the deviation around the mean value and 
vice versa. Higher CV means more resolved structures at a particular model grid 
spacing (A,) for the respective areas. A steeper slope of the curve means a more 
considerable variability of the parameter across the range of model grid spacings 
(Ar). All in all, it gives the variability of a parameter in terms of the resolved (or 
unresolved) physical structures at the given model grid spacing (Aj). 

Figure 6.2 illustrates the CV of orography (a), soil type (b) and TAI (c) of Al, 
A2 and A3, respectively. The inter-variability of resolved orography (Figure 6.2a) 
among the areas is lowest for Al (22296) and highest for A3 (33%). That means 
A3 contains more disperse orographic structures than Al. The sequential arrange- 
ment of the mean CV curves for the given set of model grid spacings of the three 
areas (Al, A2 and A3) follows the orography hierarchies arranged in the order 
of ascending mean terrain height of the regions. Here, A3 has a more significant 
slope of 15.8° than that of A2 (slope of 11.3°) and Al (slope of 3.4°) which means 
that the variability across the model grid spacings (from Asooom to A156m) for oro- 
graphy is largest for A3 and lowest for Al. Thus most of the orographic features 
are already resolved at Asooom in Al whereas Asoo0m in A3 contains a relatively 
high number of unresolved features which are yet to resolve at finer grid spacings. 
Furthermore, Al at 6125072 A625m, A312m, and Aj56m grid spacings show higher 
values of mean CV for soil type (Figure 6.2b) in comparison to A2 and A3, which 
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Figure 6.2.: Classification of the three areas and different model grid spacings based on 
areal mean Coefficient of Variation (CV) of orography (a), soil type (b), and 
TAI (c). 


means that the soil type heterogeneity is more prominent over the flat terrain (A1) 
than the other two areas (A2 and A3) except at the scales of Asooom and Aosoom 
where A2 shows slightly higher CV in comparison to Al. This is because A2 
at Asooom has a relatively larger deviation (0.4 %) of soil type (54 % sand, 0.85 
% sandy loam, 43.7 % loam, 0.85 % loamy clay and 0.71 % peat) around the 
mean value whereas Asooom in Al is dominated with similar soil type varieties 
(11 % sand, 41.2 96 sandy loam, 42.8 % loam, clay 0.4 % and 4.7 % peat) but 
smaller deviation. Finally, Al has the largest slope of 14.5" compared to A2 
(1.7?) and A3 (5.3?) which means that the soil type based variability in terms of 
unresolved/resolved soil types within a scale range from Asooom to A156m is higher 
for the flat terrain (A1) compared to that over A2 and A3 region. As such, small 
lakes and the areas with soil type peat are not resolved at Asog9m over Al, but only 
at the finer grid spacings (Aj259m to A156,,). 


40 


6.2. Convective precipitation in control runs and sensitivity experiments 


Similarly, the variability associated with TAI (Figure 6.2c) at Aj56 is highest for 
A3 ( ~ 13.5%) compared to that of Al (11.9%) and A2 (11.8%) which means 
A156m 1n A3 resolves a higher number of heterogeneous vegetated surfaces. A3 
has the highest slope of 53.6°, hence the highest difference in variability, whereas 
A2 and Al with the slope of 47.9°, 41.3° show a relatively lower level of disperse 
vegetation resolved across the given model grid spacings. Despite the relative 
differences, the range of variability of TAI across the model grid spacings remains 
significantly large for all the areas. 

From the above classification, it would be interesting to see whether this 
substantial modification of TAI values across the model grid spacing has a con- 
siderable impact on the spatio-temporal behaviour of precipitation, in comparison 
to orography and soil type. Moreover, as the selected days (Chapter 3) are char- 
acterised by locally induced precipitation, the triggering may vary significantly 
depending upon the resolved land-surface properties and the orography at the dif- 
ferent grid spacings. This dependence is discussed in detail in Chapter 7. 


6.2. Convective precipitation in control runs and sensitivity 
experiments 


This section investigates the sensitivity of spatial patterns and the areal mean of 
precipitation to different model grid spacing and land-surface resolution. 


Spatial pattern of daily accumulated precipitation in control runs 


Figure 6.3 and 6.4 show the horizontal snapshots of total accumulated precipit- 
ation until 2000 UTC from the control runs of the cases simulated over Al (flat 
terrain), dated 26 July 2012 and A3 (complex terrain), dated 12 Aug 2015, re- 
spectively (details of numerical simulations and used physics are described in 
Chapter 5). In Figure 6.3, the reference run (61567 of the control run) produces 
clustered, and non-uniform precipitation patterns in the east and southeast of Ber- 
lin whereas A1250m and the other LES grid spacings form dispersed, but relatively 
uniform-sized precipitating cells throughout the East of Berlin. Similarly, in Fig- 
ure 6.4, the reference run results into clustered precipitation patterns in the North- 
ern Black Forest along with a few small-sized precipitating cells in the south- 
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western part of the Black Forest, while, A1250m produces many much bigger and 
almost uniform sized precipitating cells in most of the area. 

Considering the comparison across model grid spacings (A,) with the same phys- 
ics, Figure 6.3c-f 1s one example showing the different behaviour of precipitation 
patterns across LES grid spacings simulated over Al, case: 26 July 2012. Ajasom 
produces intense precipitating cells in comparison to other finer LES model grid 
spacings which is a consistent model behaviour found in the precipitation patterns 
simulated for other cases (see Section 3 for details of cases) except in A2 (17 Sept 
2017) and A3 (29 May 2017) where 41567 produces stronger precipitation cells 
than other coarser LES model grid spacings (A,) (Figure B.5, B.6, B.7, B.8). 
Furthermore, the total precipitation patterns differ considerably between NWP 
(6500072 A2500m) and LES set of model grid spacings (6125077 A625m, A312m and 
A156m) in both of the cases (Figure 6.3 and 6.4). In Al, NWP model grid spacings 
produce completely different, bigger and intenser precipitating cells in compar- 
ison to LES model grid spacings (Figure 6.3) while for the case simulated over 
the complex terrain, located in the Black Forest (A3), the NWP set does not pro- 
duce precipitation at all (Figure 6.4). 
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Figure 6.4.: Same as Figure 6.3 but over A3, case: 12 Aug 2015 
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Figure 6.3.: Total accumulated precipitation patterns (colour coded) across the six model 
grid spacings (As000m, A2500m> A1250m» A625m, A312m and A156m) at 2000 UTC 
over Al, case: 26 July 2012 


Thus, a different treatment of turbulence closure can even lead to a no-precipitation 
situation. The restriction to vertical turbulent exchange and an assumption of the 
negligible contributions of a horizontal component leads to different forcing of 
convection. This assumption is strongly violated in complex terrain where 3D 
effects are important, for example, the advection and the horizontal shear pro- 
duction, act as source terms of TKE (Arnold et al., 2012), make the turbulence 
parameterization even a more sensitive parameter to choose (Honnert and Mas- 
son, 2014; Couvreux et al., 2016; Rotach et al., 2017). 

As obvious from Figures 6.3 and 6.4, there are considerable differences between 
the precipitation patterns at the transition from grid spacing Aj259 to A2500m. Al- 
though focusing on the sensitivity of convective precipitation to model grid spa- 
cing and land-surface resolution, it is worthwhile to know the impact of the differ- 
ent turbulence parameterization schemes (here, NWP and LES) on the precipita- 
tion. Therefore, to segregate the impact of model grid spacing and the transition 
of turbulence parameterizations sensitivity simulations at Aj250m for the case over 
Al (26 July 2012) are performed using the same model setup. However, applying 
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the two turbulence parameterization schemes. Figure 6.5 shows the spatial pat- 
terns of precipitation at 6412502 model grid spacing under two different turbulence 
parameterizations. A1559,, -NWP produces intenser and bigger precipitating cells 
than the smaller and scattered cells of A1250m -LES. The results of these sensitiv- 
ity simulations indicate that the differences between the precipitation patterns in 
A1250m to Azsoom (Figure 6.3) are a combined effect of the change of grid spacing 
and the turbulence parameterization alike. In the following, therefore, focus more 
on the comparison between the simulations from 4156౫ to A1250m- 
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Figure 6.5.: Total accumulated precipitation (colour coded) for Aj250m in NWP (a) and 
LES (b) turbulence parameterization schemes at 2000 UTC over Al, case: 26 
July 2012. 


Diurnal cycle of areal mean accumulated precipitation amount in 
control runs 


The diurnal cycle of areal mean precipitation calculated in the six model grid spa- 
cings is investigated. Figure 6.6 is a typical example of the time evolution of the 
area-averaged accumulated precipitation for a case simulated over Al (26 July 
2012). The precipitation in the reference run starts at 1410 UTC and ends at 1820 
UTC resulting in a total of 0.82 mm precipitation. However, the onset of precip- 
itation in Aj259m is already at 1130 UTC and it ends at 1630 UTC with the mean 
precipitation of 2.26 mm. That means in A1250m the onset of precipitation is nearly 
two hours earlier and the duration of precipitation is one hour longer than in A156,,. 
Across the LES model grid spacings, the accumulated mean precipitation amount 
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increases with increasing model grid spacing. Including the diurnal cycles of the 
precipitation of the other cases and areas (Figure B.9, B.10, B.11, B.12, B.13), a 
consistent model behaviour exists which says that the onset of areal mean precip- 
itation is earlier in Aj259 by ~1-2 hours in comparison to A156. It is also shown 
in Figure 6.6 that the mean precipitation in NWP (Aso00m) is being delayed by 
almost 1.5 hours as compared to LES (415677) model grid spacing. This is found 
as a consistent behaviour in other cases, where NWP (Aso00m) precipitation occur- 
rence is by 1-2 hours later (or even does not occur at all, Figure 6.4) than in LES 
(Ais6m) (Figure B.9, B.10, B.11, B.12, 3.13). 

The difference in onset of precipitation between the two turbulence closures 
is in agreement with the recent simulations done for A,, = 25 m to 1000 m in 
Chow et al. (2019). 
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Figure 6.6.: Time-series of total accumulated precipitation for NWP (Asogom, A2500m); 
dashed-lines, and LES model grid spacings(Aj259m, A625m, A312m, A156m)s 
solid-lines, over Al, case: 26 July 2012 


Spatial pattern of daily accumulated precipitation in land-surface 
sensitivity runs 


To address the impact of land-surface modification on convective precipitation, 
the numerical simulations have been performed by using the coarser resolution of 
land-surface properties altogether (LO_5000 m, LO_1250 m) but refined model 
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grid spacing from Asgo9m to A156m for all the six cases. See Section 5.5 for de- 
tails of the simulations strategy. Additionally, the LL. 5000 m, L. 1250 m, O. 5000 
m, O 1250 m runs have been performed only for the cases over A2 and A3 to 
understand the relatively dominant sensitivity of convective precipitation to the 
modification of land surface properties (L) and orography (O). Figure 6.7 is an 
example showing the composites of accumulated precipitation at 2000 UTC from 
the land-surface sensitivity simulations performed over A1 on 26 July 2012. The 
spatial distribution of the precipitation patterns is similar in the three runs except 
for the clustered precipitating cell in the east of Berlin, which is relatively bigger 
in size than the one in the reference run. Similarly, the other simulations with 
modified land-surface (L) and orography (O) separately from 156 m to 1250 m 
(L. 1250 m, O 1250 m) and 5000 m (L. 5000 m, O 5000 m) result into similar 
precipitation patterns as in the reference runs (Figure B.14, B.15). 


Ref A156m LO 1250m A156m LO 5000m 


(b) wp | s Æ (c) 


52.8°N} 


134E 14.2°E 13.4% 


Figure 6.7.: Total accumulated precipitation from reference run (A156,4- control run) (a), 
A156m- LO. 1250 m run (b) and Aj56,,-LO_5000 m run (c) at 2000 UTC over 
Al, case: 26 July 2012. 


Diurnal cycle of areal mean accumulated precipitation amount in 
land-surface sensitivity runs 


The impact of the coarsening of land-surface resolution on the diurnal cycle of 
precipitation is given in Figure 6.8, which shows the time-evolution of areal mean 
total precipitation calculated in the above sensitivity simulations for the same case 
over Al (Figure 6.7). The modulation in land-surface resolution over the scale 
of 1250 m and 5000 m has almost negligible impact on the onset and duration 
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of precipitation. However, the accumulated mean precipitation amount, in this 
case, is increasing with coarsening of land-surface resolution. The findings that 
the onset and duration of mean precipitation is unaffected by the modification of 
land-surface resolution also holds for the other cases (Figure B.16, B.17). 
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Figure 6.8.: Time-series of total accumulated precipitation from reference run (red square 
markers), Aj56m- LO_1250 m run (green coloured and smaller Ð) and A156m- 
LO_5000 m run (blue coloured and bigger Ð) at 2000 UTC over Al, case: 26 
July 2012. 


Area-averaged daily accumulated precipitation amount 


The overview of the precipitation derived from the simulations with different 
model grid spacing for all the cases 1s shown in Figure 6.9a. It shows the Rel- 
ative Percentage Difference (RPD) of areal mean daily accumulated precipitation 
with respect to the corresponding reference run (Aj56,,- control run of the given 
case) for the set of LES model grid spacings (A;). The areal mean daily accu- 
mulated precipitation amount of the reference runs of the different cases is listed 
in Table 6.1. In Figure 6.9a, the RPD of areal mean daily accumulated precip- 
itation is consistently decreasing with finer model grid spacing (4125077 ^e25;m. 
Aa312m; 156m) except in A2-case: 17 Sept 2017 and A3-case: 29 May 2017 where 
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A^156,, runs produce slightly higher or similar precipitation as in the other grid 
spacings. This is mainly caused by more intense precipitating cells in A156,, com- 
pared to A312m (Figure B.5, B.8). Positive (negative) RPD values mean that the 
mean precipitation at the given grid spacing is above (below) the reference run 
of the individual cases. The variability of the areal mean precipitation in 612507 
relative to the reference run is in the range of -26 to 400 % with the third quartile 
(75th percentile), Q3, of 155 % (Figure 6.92). 

Furthermore, the mean precipitation amount caused by land-surface resolution is 
smoothed down from 156 m to 1250 m results in the RPD ranging from -17% to 
37 % with the Q3 of = 7% while the range increases to -17 to 49 % and Q3 to 
22 %) when the land-surface resolution is smoothed down from 156 m to 5000 
m. Therefore, depending upon the cases, the change in model grid spacing (from 
A1250m to A156m) leads to a larger variability (~ 155%) in the areal mean precipit- 
ation amount than that caused by the land-surface modification (~ 7-22%). This 


Daily accumulated mean 


Areas Cases precipitation from reference run 
(^156,,-control run) (in mm) 
A1 26 July 2012 0.82 
17 Sept 2017 1.20 
09 June 2018 0.91 
A2 
17 Sept 2017 3.08 
12 Aug 2015 0.23 
A3 
29 May 2017 2.03 


Table 6.1.: The areal mean daily accumulated precipitation amount in the reference runs 
of different cases. 


offset in model behaviour gives a hint that there is a lack of convergence in sim- 
ulated mean precipitation or even the cloud system properties at least for model 
grid spacing >Aj56m. Bryan et al. (2003) has looked into the appropriate model 
resolution necessary for the simulation of deep moist convection. They analysed 
the squall lines simulated with grid spacings between 1 km and 125 m. The results 
show that there are no systematic trends in the fields like vertical velocity and rain- 
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water as the model resolution is increasing, and also the convergence behaviour 
of simulated squall line is still missing between 250- and 125-m grid spacings. 
This section answers the first key question of the thesis, which investigates the 
relative impact of model grid spacings and land-surface resolution on convective 
precipitation. The results reveal that the onset of areal mean precipitation is later in 
A156, than that in Aj259 (Figure 6.6). Also, in LES model grid spacings (412507 
to A156m) areal mean daily accumulated precipitation amount is highly sensitive 
to the model grid spacing and shows a stronger variability (22155 %) (Figure 6.9a) 
in comparison with the modifications of land-surface resolution, where the variab- 
ility 1s significantly smaller (227-22 96) (Figure 6.9b). This finding raises further 
questions: 


I. Why the onset and areal mean precipitation amount differ across LES model 
grid spacings? What are the processes explaining these differences? 


II. What are the most relevant factors resulting in the different areal mean pre- 
cipitation amount? 


III. What causes the varying precipitation patterns in LES grid spacings? 


The above questions I and II are elaborately discussed in the Section 6.3. To un- 
derstand the varying precipitation patterns (Question III), the resolved triggering 
mechanisms by the given LES control runs are investigated in detail in Chapter 7. 
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Figure 6.9.: RPD of areal mean daily total accumulated precipitation amount across LES 
model grid spacings (Aj259m, A625m; A312m, A1s6m) Of the control runs (LJ) (a) 
and in 41567 of LO 1250 m run (@), L 1250 m run (+), O 1250 m run (@), 
LO_5000 m run (అ. L_5000 m run (+), and O_5000 m run (o) (b) for the 
six simulated cases. Q1 and Q3 represent the first (25th) and third quartiles 
(75th percentiles). 
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6.3. Possible processes which lead to differences in precipitation: 
Clouds 


This section describes the possible processes resulting in an overall variability of 
x155 % in the mean precipitation amount between the two model grid spacings in 
LES: A1250m and A156m. The processes are discussed for the case simulated over 
Al, dated 26 July 2012, as this day has a considerable variability in the order of 
x175 % in the areal mean precipitation amount between A1559,4- and A156, COn- 
trol runs (Figure 6.9a). Figures 6.10a and 6.10b show the time-height cross section 
of area-averaged cloud liquid water content and the time series of areal mean pre- 
cipitation amount in Ajs6m and Aj259m, respectively. The spatially-averaged ver- 
tical velocity variance (1/2) represents the resolved turbulence activity at a given 
model grid spacing (A,). The Lifting Condensation Level (LCL) is based on the 
difference between the 2-m temperature and dew-point temperature, and applying 
the formula from Iribarne and Godson (2012). 

In 615677 the resolved vertical velocity variance (turbulence) (w?) within the CBL 
starts from Æ 0900 UTC followed by the first cloud formation in the CBL at Æ 
1000 UTC and the precipitation occurrence at 1410 UTC. Moreover, in À 1250m 
the resolved turbulence 1s detected at 1110 UTC, the first cloud formation at 1000 
UTC, and the precipitation at 1140 UTC. That means the resolved CBL con- 
vection is earlier in A156,, in comparison to Aj259m. Although, both of the grid 
spacings generate CBL clouds (Figures 6.10a and 6.10b) which trigger at the sim- 
ilar time, the mean precipitation is delayed by 2 hours in Aj56,, in comparison to 
A1250m- 

In the steps of investigating the questions (I and II) of the Section 6.2, the diurnal 
cycle of the cloud evolution has been analysed during its three main phases. Based 
on the cloud evolution in A156, (Figure 6.10a) the classified phases are: 


a. Environment of CI / CBL processes 
b. Cloud evolution 


c. Precipitation 
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a. Environment of CI / CBL processes 


A156m indicates that considerable resolved turbulence (w2) in the CBL already 
exists at 1000 UTC (Figure 6.11). This has an impact on the turbulent transport of 
heat and moisture in the CBL. For example, at 1000 UTC the turbulent transport 
of heat in the lowest 100 m, 1.e. in the surface layer, is mainly handled by the 
parametrized (subgrid-scale) turbulent sensible heat flux (Figure 6.12). However, 
above the surface layer the turbulent transport of heat is taken over by the resolved 
(grid-scale) sensible heat flux. The formulation for the resolved (grid-scale) and 
parametrized (subgrid-scale) fluxes, respectively are given in Appendix C.2. The 
effective mixing by resolved convective cells in the CBL results in the develop- 
ment of typical CBL profiles of 0, and q, (Figure 6.132 and b). That means the 
cloud-topped CBL consists of an unstable surface layer and a well-mixed layer 
above. 

In Aj250m, throughout the CBL there is no considerable resolved convection simu- 
lated at 1000 UTC (Figure 6.11), so that the turbulent mixing of heat and moisture 
must be done by the parametrized turbulent heat fluxes. This is obvious in the 
sensible heat transport in A1250m (Figure 6.12). The contribution by the resolved 
sensible heat flux is minimal. That means the partitioning of the total turbulent 
fluxes into resolved and parametrized flux contributions is completely different in 
A1s6m and 6125072. This also has a strong impact on the resulting profiles of 0, 
and q, (Figure 6.13a and b). On the one hand, the near-surface temperature, as 
well as the humidity, are higher in A1250m compared to Aj56,,. On the other hand, 
the resolved convection in A,56,, results in deeper and more homogeneous mix- 
ing than in Aj259 with parametrized turbulence (Figure 6.13b), where the mixed 
layer is still slightly unstable stratified at 1000 UTC (Figure 6.13a). These dif- 
ferences also have considerable consequences for the pre-convective conditions in 
the two grid spacings, as demonstrated later in this section. To understand the 
reasons for the different evolution of the environments for CI and intensification, 
the areal mean budget terms of heat (0,) and moisture (gy) are calculated. The 
budget terms evaluate the contribution of each component (for example heat ad- 
vection, phase change, radiative cooling etc.) in the prognostic equations such as 
the conservation of heat (0,) and moisture (qy). The detailed budget components 
and their formulation are discussed in Appendix C. Figure 6.14 shows the vertical 
profile of the different components of area-averaged heat and moisture budget in 


22 


6.3. Possible processes which lead to differences in precipitation: Clouds 


g kg™* 
[a] T= 
" 0.10 
1 ic 0.09 
op 0.08 
] 0.07 
= 64 0.06 
= j 0.05 
o ఖ్‌ 0.04 
i | 0.03 
| 0.02 
2+ a mum. 0.01 
J —— 7,7. 
0 EE 
€ 24- i F Ek i 
E ^, —— tot prec 
9 TA E న. 
2 54 | | | | | | | | | | | 
ee 8 09 10 11 12 Bu 15 16 17 18 19 20 
Time (UTC) 
g kg"! 
10- 
[p] 19 zus: | V 7 0.10 
I ic | 0.09 
a We EN ^ 0.08 
1 D.3 0.07 
f 6- Vv) <=] A 0.06 
= j NV మ్య 0.05 
o 42 3l » o 0.04 
m a | | 0.03 
] | LO 0.02 
2- A — | 0.01 
7 N TN 
] i s 0.00 
B pet CAS LO nt 
f= 24- i i di 09 02089899 9999797971979 979991 
E ^, —— tot prec 
9 1 Sk REINER NEM NNNM IET. cn EE AR ERRARE. | ర. ERES 
G å | | | ! | | | | | | | 
9 08 09 10 11 1? 13 14 15 16 17 18 19 20 


Time (UTC) 


Figure 6.10.: Time-height cross section of cloud liquid water content (colour coded) in 
A156m (a) and A1250m (b) -control runs for case over Al, dated 26 July 2012. 
The red and orange curves are areal mean CBL height (zi) and LCL, re- 
spectively. The black isolines indicate the spatially-averaged vertical velo- 
city variance (w/2). The black solid line with diamonds in the bottom row 
indicates accumulated area-averaged precipitation (tot prec). 
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Figure 6.11.: 
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Figure 6.12.: 
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CBL profiles of spatially-averaged vertical velocity variance (w) at 1000 
UTC in 41567 (red) and Aj259m (green) of case over Al, dated 26 July 2012. 
The horizontal lines with square and circle markers denote LCL and areal 
mean CBL height (zi), respectively. 
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Vertical profiles of areal mean mass-weighted grid-scale (solid line) and 
subgrid-scale sensible heat fluxes (dashed line) at 1000 UTC in Aj56,, (red) 
and A1250m (green) for case simulated over Al, dated 26 July 2012. The 
horizontal lines with circle and square markers denote the areal mean CBL 
height (zi) and LCL, respectively. 


1000 UTC for A156,,. The SU Me, (SUM) represents the net heating 


(moistening) rate. The MICRa, (MICR;,) represents the contribution of micro- 


physics to the net temperature (moisture) tendency, in other words, it is the net 
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Figure 6.13.: Vertical profiles of areal mean virtual potential temperature (0,) (a) and spe- 
cific humidity (qy) (b) in Aj259m (green curve) and A156, (red curve) at 1000 
UTC for case simulated over Al, dated 26 July 2012. The horizontal lines 
with circle and square markers denote the areal mean CBL height (zi) and 
LCL, respectively. 


body source terms associated with the phase changes. The TURBg, (TURB,,), 
represents the divergence of turbulent heat (moisture) flux, which describes the 
contribution of the parametrized turbulence to the net heating (moistening) rate. 
The ADVg, (ADV,,) shows the advection of temperature (moisture) by the mean 
wind. The other component is the net radiation RADg, (shown here as the sum 
of short-wave cooling, RADSWs,, and long-wave warming, RADLW6,) tendency 
terms. 

In Aj56m at 1000 UTC, the profile of the net temperature tendency, SU Me, , 
can be separated into two distinct layers (Figure 6.14a). A nearly constant net 
heating of 1 Kh“! can be found in the lowest 450 m while net cooling of up to 0.5 
Kh“! occurs between approximately 600 and 1000 m. This SU Me, profile in the 
CBL mainly results from three heat budget components: turbulent flux divergence, 
TURBg,, radiation flux divergence, RADg,, and advection of heat, ADV. 

The TURBg, contributes to the net heating in the lowest 250 m only, reaching up 
to 4.5 Kh“! close to the ground. The contribution by RAD», to net heating is 
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mainly confined to the lowest 200 m. The ADVg, contributes to the net heating in 
three different ways: cooling in the lowest 250 m, warming between 250 m and 
600 m, and cooling between 600 m and 1000 m. 

Concerning the net moisture budget, two layers can be distinguished, a layer with 
net drying of up to -0.5 gkg™! h! below 500 m, while moistening of up to 2.3 
gkg | h-! occurs between 500 m and 1000 m (Figure 6.14b). Two processes con- 
tribute to this net moisture tendency, turbulent transport of moisture, TURB,,, 
causes moistening in the lowest 250 m while ADV,, leads to drying (up to 3 
gkg 17 1) below 500 m and the moistening above. 

It is clearly seen that TURBe, (TURB,,) and ADVe, (ADV,,) are the two dom- 
inant components contributing to the SUM», (SUM,,). Further, the role of the 
areal mean ADV, term is investigated in more detail for both the temperature and 
moisture budget. ADV, is the sum of horizontal (HADV6,) and vertical advection 
(ZADVa,). ZADVa, can be decomposed into an areal mean flux and mean grid- 
scale turbulent (resolved) flux (for details of their formulations see Appendix C). 
To determine the contribution of this areal mean grid-scale sensible heat flux con- 


/ 
vergence (౨8౯) to 0164006. the amount was calculated using the model output 


fields. Figure 6.15 shows the areal mean heat advection, ADV, , and the con- 
vergence of the vertical grid-scale sensible heat flux for A156,, at 1000 UTC. As 
evident, the convergence of the vertical grid-scale sensible heat flux corresponds 
mainly to the profile of the mean advection of heat. The same finding holds for 
the convergence of the vertical grid scale latent heat flux and mean advection of 
moisture (not shown). That means the resolved convective cells transport heat and 
moisture from the lower layers (below 250 m) to upper layers (Figure 6.14a,b), 
resulting in the well-mixed layer shown in Figure 6.13a and b. Two other things 
are worth to be mentioned, too. First, the layer with negative net temperature 
tendency and positive net moisture tendency between 600 and 1000 m indicates 
the growth section of the CBL (Figure 6.14a,b). Second, as a consequence of this 
CBL growth, the areal mean CBL top reaches the areal mean LCL at around 1000 
UTC so that the CBL clouds could develop in the investigation area. 

The vertical profiles of area-averaged heat and moisture budget terms in 612507 at 
1000 UTC are shown in Figure 6.16. The overall behaviour of the net temperature 
and moisture tendency is similar to that in Ajs56,,. The contribution by radiative 
heating to the net temperature tendency is again mainly restricted to the lowest 
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Figure 6.14.: Vertical profile of areal mean heat (a) and moisture (b) budget terms within 
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Figure 6.15.: 


the CBL in Aj56,, at 1000 UTC for case over Al, dated 26 July 2012. The 
horizontal line in orange and cyan color shows the LCL and mean CBL 
height (zi), respectively. 
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Vertical profile of areal mean total heat advection (solid) and grid scale ver- 
tical sensible heat flux convergence (dashed) at 1000 UTC in Aj56,, of case 
over Al, dated 26 July 2012. The horizontal lines with circle and square 
markers denote the areal mean CBL height (zi) and LCL, respectively. 


200 m. In contrast to A156, however, TU RBe, and TURB,, dominate the temper- 


ature and moisture tendencies, respectively. That means parametrized turbulence 
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is responsible for mixing processes in the developing CBL and grid-scale flux 
convergence is of minor relevance. This is also documented in Figure 6.17, where 


w Oe and ADV», are compared. Finally, in Aj259 the mean CBL depth exceeds 
the areal mean LCL so that boundary layer clouds could develop at this time. The 
CBL clouds, which occur in A156,, and Aj259m at around 1000 UTC developed 
differently in the following hours (Figure 6.10a and 6.10b). The reasons for that 
are investigated in the next section. 

As outlined in the previous discussion, the partitioning of total turbulent heat 
and moisture fluxes into resolved and parameterized fluxes results in warmer and 
moister conditions in the lowest 500 m in Aj259 than Ajs6m (Figure 6.13). This 
affects the potential energy of the air parcels. Figure 6.18 shows the time evolution 
of area-averaged Convective Available Potential Energy of mean surface layer par- 
cel (MLCAPE) and Convective Inhibition of mean surface layer parcel (MLCIN) 
for the same case as before. The MLCAPE and MLCIN use the thermodynamical 
properties (mean temperature and humidity) of the representative shallow surface 
layer (SO hPa Æ 500 m). Aj259 calculates the maximum areal mean MLCAPE 
in the order of 221350 J kg"! at 1010 UTC and the minimum areal mean MLCIN 
in the order of +1 Jkg~! at 1110 UTC. In A1250m the MLCAPE removal occurs 
as soon as the air parcels overcome the associated MLCIN. Furthermore, A156m 
calculates the maximum areal mean MLCAPE of 21000 J kg"! at 0950 UTC and 
the minimum MLCIN in the order of +2 J kg"! at 1300 UTC. 

Higher MLCAPE (#1350 J kg!) in Ai250m together with similar MLCIN (#1 
J kg!) anticipates a higher potential to trigger deep convective systems than in 
A156m With relatively lower MLCAPE (+1000 J kg^!). This could be one factor 
explaining higher amount of precipitation in Aj259,, than in Aj56,, (Emanuel, 
1994). 

Interestingly, the MLCAPE removal in A156m already starts at 0950 UTC while 
the MLCIN is yet not overcome and before the precipitation could cause a reduc- 
tion of MLCAPE (first precipitation occurs later at 1410 UTC). The MLCAPE 
removal can be explained by analysing the MLCAPE behaviour at two times, (1) 
when peaks to the maximum value 1000 J kg! (at 0950 UTC) and (ii) at 1200 
UTC when it is reduced down to 705 Jkg"!. Figure 6.19 shows the vertical 
profiles of specific humidity (g,) within the CBL at 0950 UTC and 1200 UTC. 
At 1200 UTC the lowest ~500 m of the CBL is drier than at 0950 UTC which 
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Figure 6.16.: Same as Figure 6.14 but in A1250m 
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Figure 6.17.: Same as Figure 6.15 but in A1250m 


attributes to a drop in dew-point temperature. This breaches the resultant Level of 
Free Convection (LFC) and the rise in the LFC results into a lower MLCAPE at 
1200 UTC. That means the strong the CBL growth with deeper mixing of mois- 
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ture between 0950 UTC and 1200 UTC is the major factor for the removal of 
MLCAPE already before the precipitation occurrence. 
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Figure 6.18.: Time-series of areal mean MLCAPE for Aj259 (solid green) and Aj56,, 
(solid red) and MLCIN for A1250m (dashed green) and Aj56,, (dashed red) 
for case simulated over Al, dated 26 July 2012. The vertical dashed and 
solid lines show the time of formation of first cloud and the first occurrence 
of precipitation in the domain, respectively. 
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Figure 6.19.: Vertical profiles of areal mean specific humidity in 41567 at 0950 UTC (red) 
and 1200 UTC (blue) for case simulated over Al, dated 26 July 2012. The 
horizontal lines with circle and square markers denote the areal mean CBL 
height (zi) and LCL, respectively. 
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b. Cloud evolution 


A156m and Aj259m both result in cloud-topped CBL (Figure 6.10). However, e.g. 
until 1200 UTC clouds grow deeper in Aj250m (8000 m) in comparison to A156;; 
(226100 m). This raises the question, which processes are responsible for this dif- 
ferences in cloud development. One reason could be the higher potential energy 
in terms of favourably higher MLCAPE and lower MLCIN values in Aj259,, than 
that in Aj56, as discussed in the previous section. Additional insights could be 
gained by investigating the cloud properties and heat and moisture budgets for 
the period when CBL clouds develop into deep convection. The diurnal cycle of 
the column integrated cloud liquid water content (tqc) and water vapour (tqv) is 
shown in Figure 6.20. A1250m has the maxima of tqc in the order of 220.25 kg m~? 
at 1210 UTC in conjunction with a reduction in tqv which means that the most of 
the water vapour (g,) at this hour has been converted to cloud liquid water (qe). 
This can be also seen with the cross-section of qc at around 1210 UTC. in the 
Figure 6.100. Aj56 simulates the maxima of tqc in the order of ~0.17 kgm? 
between 1400 and 1500 UTC with only a slightly reduce increase of the mean tqv. 
That means only some part of the water vapour q, has been transformed in cloud 
liquid water (gc). With this, Aj259,, calculates higher tqc than A156,, at around 
1200 UTC which is also evident in the time-height cross section of cloud liquid 
water (Figure 6.10a and 6.10b). Figures 6.21 and 6.23 show the vertical profiles of 
the components of area-averaged heat (a) and moisture (b) budget in the CBL and 
free troposphere at 1200 UTC for Aj259 and 615677 respectively. This time was 
selected because A156,, shows no precipitation whereas Aj259,, simulates some 
precipitation (Figure 6.10a and 6.10b). This makes it a suitable hour for analys- 
ing the processes causing differences in the cloud evolution and the precipitation 
onset. 

In Aj250m, the SU Me, profile can be separated mainly into four layers: cooling 
in the lowest 100 m, warming between 100 m and cloud base, which corresponds 
with the LCL, warming from cloud base up to approximately 3 km (Figure 6.21a). 
In the layer between 3 km and 8 km, the SUMg, is small. To this SU Me,, the 
TU RBs, contributes only in the subcloud layer by heating, decreasing nearly lin- 
early from about +1.9 Kh”! at the surface to 0 Kh“! at the cloud base. Appre- 
ciable heating by radiative processes (RADg,) of up to 0.7 K h”! is restricted to the 
lowest 100 m. The most prominent contributions to the heat budget are advective 
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Figure 6.20.: Time-series of areal mean total column integrated cloud liquid water (tgc) 
(solid line) and total column integrated water vapour (tqv) (dashed line) in 
A^1565 (red) and Aj259m (green) for case simulated over Al, dated 26 July 
2012. The vertical dashed and solid lines show the time of formation of first 
cloud and the first occurrence of precipitation in the domain, respectively. 


and microphysical processes. In the subcloud layer, MICRg, causes cooling (ap- 
proximately -0.6 Kh~'), which is due to evaporation of rain, and heating by con- 
densation from the cloud base to cloud top at about 8 km. Condensational heating 
is highest with values of about 1.5 K^! in the lower part of the cloud layer. The 
contribution by ఉగ, to the SU Me, can be separated into three layers: cooling of 
up to -2.9 Kh! in the lowest 150 m, heating between 150 m and cloud base and 
cooling from the cloud base to the cloud top. Once again separating the advection 


/ 
term into the convergence of the grid-scale vertical sensible heat flux CHH and 


the remaining part of ADVg, leads to more insight into the active advection pro- 


cesses. The result is shown in Figure 6.22. In the CBL, the a profile nearly 
equals ADVg, profile. That means in A1250m at 1200 UTC, convection is large 


enough to be resolved so that the grid-scale vertical sensible heat fluxes transport 
heat from the surface layer to upper parts of the CBL. Note this roughly corres- 
ponds to the conditions of A156, two hours earlier (Figure 6.15). In the cloud 


/ 
layer, however, the contribution of - and to ADVs, , apart from the lower 500 m of 


the cloud layer, is small. That means horizontal advection is responsible for cool- 
ing in most parts of the cloud layer and by this 1s compensating the condensational 
heating. The different contribution by vertical and horizontal heat exchange in the 
CBL and free troposphere can be understood by this: in the CBL, the gradients of 
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n^ 
temperature are strongest in the vertical which makes ZADVs, respectively E : 


the decisive factor for the heat exchange. In the free troposphere, horizontal tem- 
perature differences could be considerable, too, as temperature differences exist 
between clouds and their environment. The moisture budget of Aj259 is shown 
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Figure 6.21.: Same as Figure 6.16 but for CBL and free troposphere and at 1200 UTC. 


in Figure 6.21b. Concerning the net moisture tendency, two main layers exist: dry- 
ing in the subcloud layer and moistening above - interrupted by a shallow layer 
with drying just on top of the CBL. Moistening by the TU RB,, is restricted to the 
lowest 100 m. In the subcloud layer, the evaporation of rain contributes to a slight 
increase of moisture (approximately +0.1 to 40.2 gkg ! A^). In contrast, in the 
cloud layer, condensation contributes to a reduction of moisture, with maximum 
values of -0.5 gkg~!h7! in the lower part of the cloud layer. Finally, ADV,, dom- 
inates the net moisture tendency profile. Similar to the heat budget, investigating 
the convergence of the grid-scale vertical latent heat transport with ADV, , shows 


that Ð corresponds to ADV,, (Figure B.18). That means grid scale vertical 
latent heat fluxes mainly transport moisture from the CBL into the lower cloud 
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Figure 6.22.: Same as Figure 6.17 but for CBL and free troposphere and at 1200 UTC. 


layer. Additionally, way moistens the upper part of the cloud layer, too, which 


is possible because the mean humidity profile shows a decrease with height. 

Unlike the 412507 the A156, results into a net heating layer in the subcloud layer, 
caused by the same set of the processes of the heat budget which were already 
active at 1000 UTC (Figure 6.14a) i.e. TU RBe,, RADe,, and ADVg,. The cooling 


by ADVs, in the lowest 200 m of the CBL and warming up to the cloud base is 


again due to the resolved convection (Mo) (Figure 6.24). The moisture budget 


is also still very similar to the one at 1000 UTC (Figure 6.14b), TURB,, causes 


moistening of the lowest 200 m of the CBL and resolved convection (252) 


transports the moisture into the upper part of the CBL (subcloud and cloud layer) 
(Figure 6.23b, B.19). Between 1000 UTC and 1200 UTC, the CBL clouds grew 
already deep into the free troposphere reaching up to ~ 6100 (Figure 6.10a). In the 
cloud layer, the budget terms are as follows: the net heating (0.4 Kh7') is limited 
until +3000 m only. This heating has a positive contribution from MICRae, which 
is partly compensated by a simultaneous advective cooling from ADVg,. There is 
simultaneous moisture consumption in this layer indicating cloud formation in the 
lower half of the cloud layer. On top of the net heating layer, there exists a net 
cooling layer (-0.5 K^!) where MICRg, shows the dominant contribution. This 
indicates the existence of evaporative cooling and it results into a simultaneous 
release of moisture in this layer (Figure 6.23b). Consequently, there are some 
major differences between the heat and moisture budget of Aj56,, and A1250m at 
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Figure 6.23.: Same as Figure 6.14 but for CBL and free troposphere and at 1200 UTC. 
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Figure 6.24.: Same as Figure 6.15 but for CBL and free troposphere and at 1200 UTC. 


1200 UTC (Figures 6.23 and 6.21). While cloud formation (M/CRg¢,) in A1250m 
reaches from the cloud base up to the cloud top in 8000 m, in A156,, MICRg, 
contributes to cloud formation in the lower half of the cloud layer (3000 m) while 


65 


6 Impact of model grid spacing and land-surface resolution on clouds and precipitation 


contributing to cloud dissolution in the upper part of the cloud layer (3000 to 5500 
m). 

To further understand the clouds behaviour in this evaporative cooling layer 
between 3000 and 5500 m, the spatial distribution of M/C Rg, at a height of ~4000 
m is shown for Aj250m and Aj56m in Figure 6.25. It shows a zoomed view of 
typical clouds for the two grid spacings. Positive MIC Rg, indicates Latent Heat 
Release (LHR) at the cloud core while negative M/CRg, is an indication of asso- 
ciated cooling. It is apparent that 412507 produces bigger clouds whereas 61567 
simulates more number of small clouds. Though the heating due to M/C Rg, in the 
cloud core is comparable in both of the grid spacings, the clouds in 61567 show 
significantly stronger cooling at the cloud edge and shell region in comparison to 
A1250m- The strong evaporative cooling in A156, between 3000 and 5500 m causes 
the dissolution of clouds so that on average the clouds do not grow that deep at 
that time (Figure 6.10a). The possible cause behind the occurrence of this cooling 
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Figure 6.25.: Spatial distribution of temperature microphysics tendency (M/CRg,) (colour 
coded) in A1250m (a) and Aj56m (b) at ~ 4000 m height at 1200 UTC for case 
over Al, dated 26 July 2012. The black contours are the isolines for cloud 
liquid water content of 0.1 gkg~! at that level. 


layer is explained using the size distribution of clouds calculated for the cloud 
cover simulated in A1250m and Aj56m. Figure 6.262 and b are snapshots showing 
vertical cross section of equivalent diameter of the simulated cloud cover and the 
numbers of clouds formed at 1200 UTC in A1250m and A156m, respectively. Equi- 
valent diameter equals to the diameter of a circle whose surface area 1s equal to the 
respective contour area. As evident, Aj56 grid spacing produces more number of 
smaller and shallow clouds than the bigger and deeper clouds of A1250m. Morton 
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et al. (1956) and Turner (1963) have given a quantitative description of entrain- 
ment which is based on a laboratory water tank experiments of thermal plumes, 
describing the relation between mass flux (M), entrainment rate (€) and the radius 
of the rising plume (R) (Eq. 6.2). 


— ed (6.2) 


It says that the entrainment rate is approximately inversely proportional to the 
radıus of the rısing plume. That means the smaller the plume, the higher the en- 
trainment rate. Therefore, the smaller clouds in 1567౫ are easily get detrained by 
the dry environmental air, which results in the dissolution of deep clouds in Ajsem 
in the early stage. The entrainment causes evaporative cooling, which further 
reduces the buoyancy of the convective parcels. The evaporation between 3000 
and 5500 m results in a moisture increase in that layer (Figure 6.23b). Moreover, 
the addition of moisture to the environment reduces the efficiency/extent of cloud 
dissolution by lateral entrainment processes. There is considerable dependence on 
the development of deep convection on the water vapour content in the middle tro- 
posphere (Metzger et al., 2014). The evaporative cooling in the small-sized clouds 
is a prominent phenomenon analysed in Ajs6m grid spacings also in the other case 
(A2, case: 09 June 2018), resulting in the formation of detached shallow-cumulus 
clouds (Figure B.20, B.21). 

In addition to the latent heating and the evaporative cooling, the mean up- 
draught and downdraught speed inside the clouds and its environment provide 
macro-physical characteristics of clouds. Figure 6.27 shows the vertical profiles 
of areal mean MICRg,, mean updraught in the cloudy regions (defined by qe > 0), 
mean downdraught speed in the cloud- edge and shell regions, and the normalised 
grid fraction of cloudy region for A1250m and A156,;, at 1200 UTC. The shell region 
is taken as the cloudy and non-cloudy regions having vertical velocity w < -0.5 
ms |. 

As discussed in Section 6.3, at 1200 UTC the clouds are precipitating in Aj250m 
but not in Ajs6m. In A156,,, on average heating due to MICRg, occurs from the 
cloud base up to approximately 3000 m and evaporative cooling from 3000 to 
5500 m (Figure 6.272). Grid fraction with clouds is highest just above the cloud 
base (5 96) decreasing nearly continuously up to the cloud top (Figure 6.27d). In 
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Figure 6.26.: Vertical distribution of equivalent diameter of clouds and number of clouds 
(colour coded) in A1250m (a) and A156m (b), respectively at 1200 UTC for 
case over Al, dated 26 July 2012. 


the cloud layer, the mean updraught speed increases from the cloud base up to 4.25 
km (Figure 6.27b), which corresponds with the level where evaporative cooling 
is strongest (Figure 6.27a). The mean downdraught speed also reaches maximum 
values in the upper part of the cloud (-1.2 ms !, Figure 6.27c) i.e. where the 
evaporative cooling is strongest Figure 6.27a). A typical feature in Aj56,, at 4 
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km, 1.6. the level with the lowest M/CRg, value, during the period when the cloud 
fraction grows, 1s that the donwdraughts occupy a greater area than the updraughts 
(Figure 6.28). 

In 125077 the cloud fraction is approximately 6 % just above the cloud base (Fig- 
ure 6.27d). A secondary maximum of 3.5 % can be found at about 4500 m. At 
about 5500 m, the cloud fraction decreases, reaching a value of about 1%. This 
part of the cloud layer with higher cloud fraction, 1.6. between the cloud base 
and 5.5 km, corresponds with heating due to LHR Figure 6.27a). Concerning 
up- and downdraught speeds, the values are rather homogeneously distributed 
over the cloud layer (Figure 6.27b,c) and the grid fraction clouds with up- and 
downdraughts at 4 km, at least during the growing phase of cloud fraction (1100 
to 1200 UTC), is quite similar (2 %) (Figure 6.28). 

To summarise the main differences between Aj56,, and in Aj259,, at 1200 UTC: 
in 61567 the profile of the downdraught speed exhibits considerable differences, 
being stronger in the layer where evaporative cooling exists (around 4 km) than in 
the layers with condensational heating, while in 6125077 the downdraught speed 
is more similar in the whole cloud layer. Additionally, in A156,, at 4 km the grid 
fraction with donwdraughts exhibits the one with updraughts while Aj259, it is 
nearly the same. That means the layer with evaporative cooling has a considerable 
impact on the macro-physical properties of the clouds. Therefore, the evaporative 
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Figure 6.27.: Vertical profile of area mean temperature microphysics tendency, MICRa, 
(a), updraught speed in cloudy region (b), downdraught speed in the cloud- 
edge and shell area (c), and grid fraction with clouds (d) in A156, (red) and 
A1250m (green) at 1200 UTC of case over Al, dated 26 July 2012. 
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Figure 6.28.: Time-series of mean updraught (solid) and downdraught (dashed) grid frac- 
tion at Æ 4000 m for cloudy region in Aj250m (green) and A56, (red) of case 
over AT, dated 26 July 2012. 


cooling results in the dissolution of clouds in this layer which inhibit the growth 
of clouds, and this explains why there is no precipitation in Ajs6m at 1200 UTC. 


c. Precipitation 


As obvious from Figure 6.10a, first precipitation occurs in A156,, at about 1350 
UTC. Therefore, the main question arises what happens between 1200 UTC and 
1400 UTC allowing the clouds to grow deeper and to precipitate. This is demon- 
strated in Figure 6.29. In the CBL and the cloud layer turbulence (w?) becomes 
stronger and reaches higher up, i.e. the 0.2 m*s~* isolines increased from about 
5 km at 1200 UTC to 8 km at 1400 UTC. Simultaneously, latent heat release in 
the lower part of the cloud layer (approximately 1.8 km — 3 km) intensified reach- 
ing values of MICRg,- 1.5 K h“! when the first precipitation occurred. Shortly 
after that, i.e. at about 1430 UTC, the evaporative cooling in the layer around 
5 km nearly diminishes. In parallel to the changing of the microphysical cloud 
processes, the macro-physical cloud conditions also changed. The clouds became 
broader and deeper. For example, at 1400 UTC clouds with an equivalent diameter 
of 228.75 km can be found reaching up to 5.5 km while few clouds tops reached 
up to the tropopause (Figure 6.31). This cloud growth between 1200 UTC (Fig- 
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ure 6.32a) and 1400 UTC (Figure 6.32b) is attributed with cloud aggregation. 
The cloud aggregation can be seen in the spatial distribution of low-level cloud 
cover at 1400 UTC (Figure 6.32b) where the aggregated clouds form mainly in 
the east and northeast of Berlin. Further, the height vs longitude diagram after 
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Figure 6.29.: Time-height cross-section of temperature microphysics tendency M/CRg, 
(colour coded) in Ajs56-control run for case over Al, dated 26 July 2012. 
The red and orange curves are areal mean CBL height and LCL, respectively. 
The black contours are area-averaged vertical wind variance (w^). 


averaging the cloud liquid water content over the latitude from 51.6? N to 52.97 
N at 1400 UTC (Figure 6.33) clearly shows that the aggregated clouds extend till 
5 to 6 km. Between 1200 to 1400 UTC, there is an increase of the areal mean 
relative humidity in the free troposphere up to 5.5 km (e.g. about 18 % at around 
4 km) (Figure 6.30). It indicates that the dissolution of clouds moistens the free 
troposphere, which makes the atmosphere less dry. In the subsequent 3 hours, 
convective precipitation occurred (clustered precipitation patterns in Figures 6.3 
and 6.32c) with evaporative cooling in the subcloud layer due to evaporation of 
precipitation (Figure 6.29). Until 1600 UTC, MLCAPE degraded considerably 
(Figure 6.18) and cloud liquid water (tqc) was removed completely from the at- 
mosphere (Figure 6.20) so that the convective activity stopped (Figure 6.29). 
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Figure 6.30.: Vertical profiles of mean relative humidity in Aj56,, at 1200 UTC (red) and 
1400 UTC (blue) for case simulated over Al, dated 26 July 2012. 
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Figure 6.31.: Vertical cross section of cloud size distribution in Aj56,, at 1400 UTC for 
case over Al, dated 26 July 2012. 
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Figure 6.32.: Low-level cloud cover at 1200 UTC (a), at 1400 UTC (b), and total accu- 
mulated precipitation until 1500 UTC (c) in A156m of case over Al, dated 26 


July 2012. 
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Figure 6.33.: Latitudinally averaged time-longitude diagram of cloud liquid water content 
at 1400 UTC in the first row and total accumulated precipitation at 1500 
UTC in the second row in A156, of case over Al, dated 26 July 2012. 


Lesson drawn 


The first take-home message here is that the areal mean daily accumulated pre- 
cipitation amount in LES model grid spacings (612507 to Aj56m) is more sensitive 
to the model grid spacing and shows a stronger variability (+155 %) in compar- 
ison to the land-surface resolution (7-22 96). Across the LES grid spacings, on 
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most of the investigated days the areal mean precipitation amount decreases with 
increasing grid spacing. Additionally, precipitation onset is delayed by 1 or 2 
hours for the higher grid spacing runs. This is associated with the fact that in the 
simulations with finer horizontal grid spacing, first smaller clouds develop which 
show intensive evaporative cooling at their shell and edges. This restricts the cloud 
growth and therefore suppresses precipitation. Later on, cloud aggregation is an 
important process to allow these clouds to grow deeper and to precipitate. 
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This chapter describes the triggering processes of convection depending on land 
surface identified in the above-mentioned case studies. The basic ingredients of 
the triggering of convection or CI are atmospheric instability, sufficient mois- 
ture supply and wind convergence (Doswell III, 1987; Chen and Avissar, 1994; 
Kalthoff et al., 2011; Barthlott and Kalthoff, 2011; Adler et al., 2011). The sur- 
face layer which in direct contact with the Earth’s surface is influenced by the 
surface sensible and latent heat fluxes. The turbulent heat flux further transports 
heat into the CBL and causes differential heating in the atmospheric layer when 
the surface is heterogeneous. As a result, the evolution of the CBL itself experi- 
ences a diurnal cycle. According to Segal and Arritt (1992), the daytime spatial 
heterogeneity in surface sensible heat flux can lead to thermally-induced circula- 
tions. The spatial distribution of surface sensible heat flux depends on different 
factors like landforms, land use variation, soil moisture heterogeneity and cloud 
cover. That means the surface sensible heat flux is one measure to quantify the 
land-surface heterogeneity. Additionally, the triggering mechanisms could be dif- 
ferent depending upon the scale of land-surface heterogeneity and the synoptic 
situation. 

Here, the land-surface heterogeneity has been broadly categorised into land- 
surface properties and orography. The first investigation area, Al, could show 
the formation of moist convection due to variation in land-surface properties, for 
example, vegetation, soil type, land use (including water, sealed surface). The 
inhomogeneities in land-surface properties cause differential surface fluxes and 
form pressure gradients which trigger mesoscale circulations. Another key factor 
in this regard is the orography. For areas A2 and A3, additional to the land-surface 
properties, orography becomes gradually more important. 

In general, the triggering mechanisms over the heterogeneous regions could be 
based on dynamically as well as thermally-driven circulations. Under fair weather 
conditions, the thermally-driven circulation over heterogeneous terrain (e.g. at 
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crestlines and adjacent land surfaces) poses a reversal of wind direction twice 
every day (Whiteman, 2000). Subsequently, the temperature difference results in 
the density variation of air masses, which induces a horizontal pressure gradient 
and initiates the circulation. The air mass close to the surface flows from areas of 
higher pressure to the areas of low pressure. Conceptually, the upper branch of the 
current circulation runs in the opposite direction. Therefore, baroclinicity tends 
to generate mesoscale wind systems. Some examples of these circulations are 
lake breezes, slope-wind systems and valley-wind systems. Triggering processes 
which are based on differential heating and the local wind circulations are called 
as primary triggering mechanism. 

Further, the above factors precondition the atmosphere for the generation of 
moist convection which then can result in precipitation. The precipitation itself 1s 
a decisive factor which controls the formation of new convective cells. The trig- 
gering of new cells depends upon the transient primary triggering mechanisms. 
Therefore such processes are called secondary triggering mechanisms. For ex- 
ample, the triggering through gust fronts and the outflow boundaries of existing 
or decayed precipitation systems could generate new convective cells. The evap- 
oration of rain causes a cooling effect in the surface layer accompanied by a strong 
downdraught. This evaporatively cooled region of downdraught air is called a cold 
pool. This cold and dense air then spreads laterally at the surfaces and tends to lift 
the surrounding warm and light air. This process 1s called mechanical lifting (Torri 
et al., 2015). The lifting can also occur due to the convergence of gust fronts of 
the cold pools and the environment wind. Due to mechanical lifting, a cold-pool 
system is very likely to generate a small convective cell at the cold pool edge. The 
details of the genesis and formation of new convective cells due to cold pools will 
be discussed in Section 7.2.1. Moreover, the high-resolution simulations allow us 
to understand the small features embedded into cold-pool evolution for example 
associated gust fronts, generation of new convective cells. 

Another way of categorising the triggering is through the underlying sources 
of heterogeneity itself. These sources of triggering can be represented under two 
categories, namely static and dynamic. Effects under static heterogeneity com- 
prise the phenomenon generated due to surface properties variation, for example, 
lake-breeze, Urban Heat Island (UHI) effect and orographically-induced circula- 
tions (Kang and Bryan, 2011; Rieck et al., 2014). In contrast the dynamic het- 
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erogeneity pertains the processes which are triggered due to the modulation from 
cloud radiative effects or cold pools (Marsham et al., 2007; Rieck et al., 2015). 


7.1. Primary triggering 


The primary triggering mechanism can be broadly divided into two categories. 
The first category involves those locally induced processes which are caused by 
inhomogeneities in land-surface properties or due to the spatial pattern of Bowen 
ratio (the ratio of surface- sensible heat flux and latent heat flux) generated by 
heterogeneous soil type, soil moisture and vegetation, for example, thermally- 
induced circulations. The second category consists of orographically induced cir- 
culations, for example, slope winds, valley winds. The above processes occur 
independently or in a way that they are superposed. 


7.1.1. Thermally-driven circulation 
Theoretical background 


Bowen ratio gives the partitioning of the available energy into surface- sensible 
heat flux and latent heat flux. The spatial pattern of Bowen ratio modulates the 
CBL state, which impacts the conditions resulting in locally induced convective 
systems, for example, lake breeze (Section 7.1.1a) and UHI (Section 7.1.1b). A 
lake breeze can be used as an example to explain the fundamentals of a thermally- 
driven circulation system. Generally, lake breeze (sea breeze) can be seen as the 
onshore/offshore winds occurring during daytime and observed over the lake/sea. 
As the name suggests, the wind blowing from the lake (land) in the daytime (night- 
time) and blowing towards land (lake) is called lake (land) breeze. It is caused 
by the surface temperature difference of the lake and the adjacent land. This 
creates a pressure gradient at lower levels in the atmosphere, which results in the 
lake-breeze near the Earth's surface (Atkinson, 1981). The definition of lake/sea- 
breeze circulation is based on Bjerknes’ circulation theorem (Bjerknes et al., 1898; 
Thorpe et al., 2003) which is an expression obtained by taking the line integral of 
Newton's second law for a closed curve in the fluid. This gives a relation between 
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the relative circulation C (a measure of the rotation of the fluid) and the enclosed 
area projected on the equatorial plane Ag. 


R (7.1) 


Dt Dt 


In equation 7.1, the term on the left hand side denotes the rate of change of relative 
circulation, the first term on the right hand side 1s called solenoidal term. The 
solenoid term originates from the intersection of isobaric and isothermal surfaces 
due to the density variation, which results in a rotational acceleration (or solenoid). 
The second term denotes the inertial effect. Q is the angular speed of the Earth's 
rotation. Using Stoke’s theorem the equation 7.1 can be expanded as: 


DC > 
X- / (Va x Vp) dÁ, (7.2) 
A 
DC 
m. = | IVallvplisin(B) da, (7.3) 
A 


where œ and p are the specific volume (à = 5) and pressure, respectively. p 
denotes the angle between the gradients of œ and p. 

According to this theorem, the rate of change of relative circulation (C) of a 
closed chain of fluid particles is determined by two effects: a) the solenoid effect, 
in which the solenoids tend to modulate the circulation. This can be quantified 
using the total number of solenoids falling into the closed curve. This equals 
to the flux of the vector product of the gradients of each quantity (X and p in 
equation 7.2) through a surface area bounded by the curve. b) the inertial effect 
which tends to decrease the relative circulation by a quantity proportional to the 
rate of change of the enclosed area (Ae), which is projected on an equatorial plane. 

For a barotropic fluid in which the density is a function only of pressure, the 
solenoidal term is zero. That means the absolute circulation following the mo- 
tion is conserved in barotropic fluid and therefore, the relative circulation for a 
closed chain of fluid particles will be modulated only by the enclosed horizontal 
area or the average latitude changes. But unlike the barotropic fluid, the density 
of the baroclinic fluid is a function of both temperature and pressure and there- 
fore possesses a distinct air mass region, hence density gradients. Therefore, the 
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expression for Bjerknes’ theorem for baroclinic flow represents the generation of 
pressure-density based solenoid circulation, which further explains the develop- 
ment of a thermally driven lake-breeze circulation. Using the ideal gas law, the 
solenoidal term in equation 7.1 can be expressed as: 


DC, 
— QRTadl 7.4 
Di $ np (7.4) 
D S 
Ca = Rain( 2) (p m (7.5) 
Dt Di 


where C, is the absolute circulation, 7; and 7, represent the mean temperature 


over the lake and land-surface respectively. During daytime, T is greater than Tj 
DC; 
Dt 
the mechanism behind a thermally-driven circulation. 


which therefore leads to 


> 0. This resultant circulation, in general, explains 


The structure, evolution and life-cycle of sea and lake breezes are thoroughly 
studied and also reviewed both using observations (Atkinson, 1981; Simpson, 
1994; Steyn, 1998) and the numerical studies (Miller et al., 2003; Antonelli and 
Rotunno, 2007; Porson et al., 2007; Crosman and Horel, 2010). The recurrent and 
ubiquitous nature of these thermally-driven systems makes them a crucial system 
or mechanism for triggering of convection. Moreover, the same theoretical con- 
sideration is basically valid for the other thermally-driven circulations caused by 
differences in land-surface properties like UHI. UHI effect directly affects the 
mankind by generating high temperature in the urban/city areas in comparison 
to the rural or countryside areas, which influences the regional weather, increase 
the pollution level and affect the well being of inhabitants. It is the result of the 
reduction in vegetated areas, thus evaporation, or varying land use, by replacing 
green areas with roads, building and other human-made structures which absorb 
the short wave radiation more and reflect less because of their low albedo, and 
therefore, cause surface temperature to rise (Oke, 1987; Golden, 2004; Nakayama 
and Fujita, 2010). 


a. Lake breeze in ICON-LEM 


Lake-breeze systems are captured in all the LES grid spacings for the case over 
Al, dated 26 July 2012. The lake breeze is one of the triggering mechanisms 
that contributed to the formation of the convective precipitation over this area. In 
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this example, the convection is developed adjacent to Scharmiitzelsee, which is a 
lake located in Brandenburg, Germany. This lake is about 10 km long and 1.5 km 
wide surrounded by agriculture areas in the west and east, and with forest in the 
north and south. This convection is triggered and reproduced by all LES resolu- 
tions, but the clouds do not precipitate in the Aj56. As25m is used to investigate 
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Figure 7.1.: Cloud cover at 1100 UTC (a) and the accumulated precipitation until 1200 
UTC (b) in Ag25 of case over Al, dated 26 July 2012. The brown contours 
encircle the areas with lakes. 


the triggering processes over flat terrain, 1.6. the process chain from initiation 
of convection via cloud formation and precipitation. In A625m, the first cloud is 
formed at 1020 UTC in this region and the precipitation occurred at 1140 UTC. 
Figure 7.1 shows the cloud cover at 1100 UTC and the resulting accumulated pre- 
cipitation until 1200 UTC. With the north-northwesterly large scale wind, these 
clouds precipitate further to the southeast of Scharmiitzelsee. 

Furthermore, the source region of the triggered cloud and its embedded con- 
vective structure is identified using the backward trajectory. For this purpose, 
the backward trajectories are calculated using LAGRANTO (for details about 
LAGRANTO see Section 5.6) over a 50 km x 60 km box i.e. 13.85? E, 52° 
N; 14.30° E, 52° N; 13.85° E, 52.5° N; 14.30° E, 52.5° N starting from the cloud 
base at 1100 UTC (Figure 7.2). The selection criteria applied to the pre-calculated 
trajectories gives a subset of the trajectory lines originating from the cloud base. 
Figure 7.2 shows in the 3-D view of first part of the process chain. The specific 
humidity is traced back along the trajectories, which provides further insights into 
the physical properties of the air parcels. It shows that most of the air parcel 
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seeded into the clouds originated from near-surface over this period. The moist 
parcels originated above the lake, travel near-surface towards the eastern area of 
Scharmützelsee and then rise vertically, into the cloud base. Now, the surface con- 
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Figure 7.2.: Backward Lagrangian trajectories using LAGRANTO in A625,m of case over 
Al, dated 26 July 2012. The black crosses denote the starting points at the 
cloud base at 1100 UTC and the red circles denote the end points of traject- 
ory integration until 0600 UTC. The colours along the trajectories indicate 
specific humidity (¢kg~!) and the colour in the xy plane denotes orography. 


dition of the area where the parcels converge and rise are investigated. Figure 7.3a 
shows the surface sensible heat flux in background colour and the trajectory points 
which are projected in space (circles) and coloured with the specific humidity both 
at 1000 UTC. The visual inspection shows on one hand that the moist air parcels 
are positioned at the transition area of surface sensible heat fluxes before they start 
to rise. Few of the relatively dry parcels, as they are aloft to ~300 m, from the 
northwest also converge over this area. On the other hand, the areas with enhanced 
surface sensible heat flux develop in the east of the bow-shaped Scharmiitzelsee. 
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That means the region with differential sensible heat flux leads to the triggering 
of convection. 
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Figure 7.3.: Surface sensible heat flux (colour shaded) and the trajectory points (circles 
with specific humidity in the colour shade) (a), wind convergence (colour 
shaded) and wind vectors (b) at 10 m at 1000 UTC, in A625m of case over Al, 
dated 26 July 2012. The brown contours encircle the areas with lakes. 


To identify the possible factors causing the high sensible heat flux, the areal mean 
Standardized Multiple Regression Coefficient (SMRC) (Siegel, 2016) is calcu- 
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lated for surface- sensible heat flux (H) and latent heat flux (E). For example the 
formulation for H 1s: 


SMRCy = b, x —, (7.6) 
OH 

where k represents the explanatory variable like Soil Moisture Index (SMI), TAI, 
orography, net radiation, and horizontal wind at 10 m. b; is the coefficient estimate 
of H, and ర్క and og are the standard deviation of the explanatory variable and the 
outcome variable (7), respectively. These variables are selected because they 
are assumed to be the most relevant parameters determining the partitioning of 
energy exchange at the Earth’s surface into the surface- sensible and latent heat 
flux, respectively. 
The SMRC estimates how much increase in the explanatory variable affects its 
relative importance or position within the group for the determination of the out- 
come variable (7, E). The high value of SMRC for an explanatory variable shows 
its more significant contribution to the determination of the outcome variable and 
vice versa. Investigating the period before the first clouds were generated, which 
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Figure 7.4.: The areal mean SMRC for SMI in black, TAI in green, orography (Orog) in 
blue, net radiation (NetRad) in pink, and horizontal wind at 10 m (Wind10m) 
in cyan with respect to H (a) and E (b), respectively for A625m of case over 
Al, dated 26 July 2012. The vertical dotted line shows the onset time of first 
cloud in this investigation region. 


is before 1020 UTC, the sensible heat flux shows the strongest positive correlation 
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of SMRCy = +0.6 to +0.7 with the horizontal wind speed (Figure 7.4a). A much 
lower, but also a positive correlation of H exists with the net radiation (SMRCy 
~ +0.17) whereas H is only weakly negatively correlated with SMI and TAI, and 
nearly no correlation exists with orography. Moreover, before cloud formation 
occurs the latent heat flux, E, has strong correlation of SMRCg œ +0.65 is with 
TAI (Figure 7.4b). A less, but still positive correlation of E exists with 10 m wind 
speed (SMRCE ~ +0.3), followed by even a weaker dependence of E on radiation 
and SMI. Finally, the dependence of E on orography is smallest. The results 
indicate that in this case, H is not highest where the net radiation, 1.6. the avail- 
able energy for transformation into the turbulent fluxes is highest. Instead, the 
near-surface wind speed — a relevant factor in determining the vertical turbulent 
diffusion - is the dominant factor to H. It suggests that the main mechanism sup- 
porting the higher sensible heat flux is through the lake-breeze circulation system, 
which is visible in Figure 7.3b. Triggering then occurs along the sea breeze front, 
where considerable wind convergence with lifting is generated. This region is also 
corresponding with the origin of most of the trajectory points (Figure 7.3a) which 
anticipates the formation of the secondary circulation system. 


b. UHI effect 


The impact of a city on triggering of convection with subsequent precipitation also 
showed up in area Al on 26 July 2012. This effect is identified over the city named 
Fürstenwalde located in the north east of Germany (Figure 7.5). The city covers 
an area of approximately 12 km”, i.e. the diameter of the city is roughly 4 km. In 
the north, the city is surrounded by more agricultural-used land while in the west, 
south and east, forested areas dominate. The cloud of interest formed at 1110 
UTC over the southern part of Fiirstenwalde and moved with the ambient north- 
northwesterly wind further to the southeast. The cloud cover at 1130 UTC and the 
resulting accumulated precipitating system at 1230 UTC in Ag25,, in Figure 7.5 a 
and b, respectively. At 1130 UTC, the diameter of the cloud was about 4 km, and 
the cloud top reached up to 3.5 km and grew to 5.5 km before the precipitation 
started. At 1230 UTC, the precipitation in the centre of the cloud reached up to 5 
mm. 

In order to identify the triggering of this cloud system the backward trajectories are 
calculated over a box of 13.85° E, 52.2° N; 14.25° E, 52.22 N; 13.85° E, 52.5° N; 
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14.25* E, 52.5° N extent. Similar to the above example of the lake-breeze system, 
Figure 7.6 shows the 3-D view of the resulting convective system. The trajectories 
start from the grid points of the cloud base, ~1000 m, at 1130 UTC and then integ- 
rated backwards until 0700 UTC. It shows that most of the air parcel seeded into 
the clouds originated from the near-surface region over this period, and converge 
over the western premises of the city (Figure 7.7b). The surface sensible heat flux 
at 1030 UTC, 1.e. 40 min before the cloud forms, together with the projected tra- 
jectory points, are shown in Figure 7.7a. Because of the sealed surfaces in the city 
area, the surface sensible heat flux reached values of up to 220 W m ^ while in 
the surrounding, the fluxes are much less (H — 100W m ?). In the south of the 
city, the northern tip of the Scharmützelsee with its characterised low sensible heat 
fluxes (H ~ —20W m ?) over the water body and enhanced sensible heat fluxes 
east of the lake are visible again. The spatial distribution of the latent heat flux is 
vice versa (not shown), i.e. evaporation is lower over the city and higher over the 
surrounding areas. At 1030 UTC, according to the trajectory points, the air parcels 
already reached the border of the city or the city centre itself. Afterwards, they are 
lifted to the cloud base (Figure 7.6). To investigate the main dependencies of the 
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Figure 7.5.: Cloud cover at 1100 UTC (a) and the precipitating cell at 1200 UTC (b) in 
Ao25m Of case over Al, dated 26 July 2012. The magenta coloured solid 
lines enclose the area of Fürstenwalde. The investigated cloud is indicated by 
crossed lines. 


surface sensible and latent heat fluxes, the corresponding SMRCs are calculated 
again (Equation 7.6). In the period before the cloud formation, 1.6. at 1110 UTC, 
the TAI and the horizontal wind speed at 10 m show considerably stronger influ- 
ence on H in comparison to the net radiation, orography and SMI (Figure 7.8a). 
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Figure 7.6.: Backward Lagrangian trajectories using LAGRANTO in A625,m of case over 
Al, dated 26 July 2012. The black crosses denote the starting points at the 
cloud base at 1130 UTC and the red circles denote the end points of traject- 
ory integration until 0700 UTC. The colours along the trajectories indicate 
specific humidity (g kg!) and the colour in the xy plane denotes orography. 


As expected H shows a considerable negative correlation (= -0.5) with TAI and 
a positive correlation (~ +0.3 to 0.4) with wind speed at 10 m while H is only 
weakly correlated with the other variables (~ +0.1). Moreover, E shows a strong 
positive correlation with TAI and 10 m wind speed in the order of — 40.7 and — 
+0.4, respectively, while E is weakly correlated with the other factors during this 
period (Figure 7.8b). That means the reduction in vegetation or plant cover (TAT) 
in the city results into higher H, i.e. more warming of surface in comparison to the 
surrounding. This differential heating creates a horizontal pressure gradient and 
the moist parcels coming from the surrounding converge over the areas of high 
sensible heat flux (Figure 7.7a) and trigger the convection. Another example of 
triggering process based on the soil type heterogeneity is identified over A3 for 
the case dated 12 Aug 2015 (Figures B.23, B.24, B.25). 
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Figure 7.7.: Surface sensible heat flux (colour shaded) and the trajectory points (circles 
with specific humidity in the colour shade) (a), wind convergence (colour 
shaded) and wind vectors (b) at 10 m (b) at 1030 UTC, in A625m of case over 
Al, dated 26 July 2012. The magenta coloured solid lines enclose the area 
of Fürstenwalde. The negative sensible heat flux ~ —20W m”? indicates the 
Scharmützelsee. 


The above section describes the identified land surface-based triggering mech- 
anisms using an example of one of the grid spacings (Ag25m) of the control runs 
(for details of the simulation setup see Section 5.5). Triggering of the two precip- 
itating cells was observed in all the LES grid spacings of control runs (A156; to 
A625m), although the cells were slightly different in shape. The question is, how- 


87 


7 Identified triggering mechanisms 


—e—Wind1l0m 


Time (UTC) Time (UTC) 
(a) (b) 
Figure 7.8.: The areal mean SMRC for SMI in black, TAI in green, orography (Orog) in 
blue, net radiation (NetRad) in pink, and horizontal wind at 10 m (Wind10m) 
in cyan with respect to H (a) and E (b), respectively for A625m of case over 


Al, dated 26 July 2012. The vertical dotted line shows the onset of first cloud 
in the investigation region. 


ever, how does the modulation in land-surface resolution impact the triggering of 
these processes. To investigate this, an overview comparison of the lake breeze 
and UHI generated clouds is provided for the land-surface sensitivity experiments 
i.e., LO 1250 m (Figure 7.9d-f) and LO. 5000 m (Figure 7.9g-1) with respect to 
control run (A625,m, Figure 7.9a-c). With respect to the sensible heat flux patterns, 
the results show that the lakes and the city are still resolved in the LO 1250 m 
simulations - at least their substantial parts (Figure 7.9d-f) - but not in LO 5000 
m run (Figure 7.9g-1). This is, because the spatial heterogeneity of these land- 
surface structures (water bodies, city) range in between 4-5 km, for which on one 
hand the land-surface resolution of 1250 m 1s sufficient, and on the other hand 
at this spatial heterogeneity scale thermally induced circulations often tend to be 
generated (Shuttleworth, 1991). However, as the land-surface resolution 5000 m 
does not resolve these patterns anymore, secondary-circulations systems are miss- 
ing in the LO 5000 m simulations as well as the corresponding precipitating cells 
(Figure 7.9g-1). Nevertheless, precipitating cells of similar size and intensity are 
generated in other areas 1.e. other triggering processes are active because the envir- 
onmental conditions are favourable in the investigation area. This, in turn, agrees 
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Figure 7.9.: Spatial patterns of surface sensible heat flux at 1000 UTC, low-level clouds 
1130 UTC, and accumulated total precipitation at until 1230 UTC in A625m of 
the control run (a-c), LO 1250 m (d-f), and LO 5000 m (g-1) of case over Al, 
dated 26 July 2012. The clouds and precipitating cells marked with horizontal 
and vertical bars denote the respective cloud and the resulting precipitation 
through the lake-breeze circulation and UHI effect, respectively. The brown 
solid lines enclose the areas with lakes. 


with the finding that the areal mean precipitation of LO 5000 does not differ that 


much from the corresponding control run (Ag25m) (Section 6.2). 
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7.1.2. Orographically-induced circulation 
Theoretical background 


Orography may generate the dynamically-driven clouds through orographic lift- 
ing, which is needed to overcome the CIN and release the CAPE. This process 
depends upon the shape and size of the obstacle (mountains) and the background 
wind. Moreover, the thermally driven daytime upslope and up-valley winds in the 
lowest atmospheric layers can occur. The genesis of these processes are well stud- 
ied and identified (Steinacker, 1984; Atkinson, 1981; Markowski and Richardson, 
2011; Zardi and Whiteman, 2013). The prerequisites are calm and cloud-free 
conditions so that the daytime absorption of radiation by the sloping surface cre- 
ates warming of air close to the surface (Whiteman, 2000). By that, the air near 
the sloping surface becomes warmer than the air in the free atmosphere at the 
corresponding height above sea level. Based on the hydrostatic principle, a ho- 
rizontal pressure gradient will be induced. This generates a low-level horizontal 
pressure gradient, under which upslope wind occurs. Further, valley winds are 
thermally-induced circulations. They are primarily driven by pressure gradients, 
induced due to along-valley temperature differences Whiteman (1990). These dif- 
ferences can be explained by the so-called geometric or Topographic Amplifica- 
tion Factor (TAF) (e.g., Steinacker (1984); McKee and O’ Neal (1989); Whiteman 
(2000)). The volume of the air in a box above a valley atmosphere is less than 
over flat terrain. That means during the daytime, the same surface turbulent heat 
and moisture fluxes would result in a greater increase of heat and moisture in the 
air over the valley than over the adjacent foreland and results in up-valley wind 
systems. 

The following section discusses typical examples of thermally-driven circulation 
systems which are orography based and superposed with others. Unlike for the flat 
terrain cases where Ag»s,, is sufficient to show the land-surface properties based 
circulations, the 641567 results are shown for the examples of orography-based 
triggering processes as it comprehensibly represents the fine orographic based 
circulations. 
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Superposition of slope- and soil-type heterogeneity-based 
circulation 


An example of such an orographically-induced circulation which trigged convec- 
tion identified in a case simulated over the orographic area A2 dated 09 June 2018 
is discussed (details of the case is discussed in Section 3.2). On this day, the 
convection formed over the Harz mountain range. The day has east-southeasterly 
large-scale wind, and the first cloud was formed at 0920 UTC. over the mountain 
ridge. Figure 7.10 shows the orography, soil type, cloud cover at 1010 UTC and 
the resulting accumulated precipitation of the system until 1200 UTC. The oro- 
graphy indicates that the highest peaks are situated in the western part of the Harz 
mountain, like its highest peak, the Brocken (1141 m) (Figure 7.10a). The dom- 
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Figure 7.10.: Orography (a), soil type (b), low-level cloud cover at 1010 UTC (c) and the 
resulting accumulated precipitation at 1200 UTC (d) in Aj56 of case over 
A2, dated 09 June 2018. 


inating soil type in the Harz is loam (Figure 7.10b). However, in the eastern part 
of the peaks is loamy clay, while a stripe of peat stretches to the west. In the flat 
areas around the Harz mountain range, sand and sandy loam is the dominating soil 
type. At 1010 UTC, convective clouds are mainly present over the highest peaks 
of the Harz mountain (Figure 7.10c). Further, Figure 7.10d shows that this is the 
region where precipitation occurred in the subsequent 2 hours. The precipitation 
pattern covers most of the area of the Harz mountain ridge. 
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The conditions, which led to the simulated precipitation are as follows: the sens- 
ible heat flux distribution at 0830 UTC shows higher values over the Harz than 
over the surrounding flat land (Figure 7.11a). Maximum values of up to 280 
Wm? can be found over the spot with peat in the west of the Brocken. Ana- 
lysis of the dependence of the turbulent fluxes on the aforementioned parameters 
indicates that H is considerably correlated with orography (7 +0.40) and SMI (~ 
-0.38), which is also reciprocated with a considerable correlation of E with oro- 
graphy (œ -0.42), SMI (œ> +0.40), and TAI (~ +0.38) (Figure 7.13). That means 
the orography together with the soil type results in the simulated surface heat flux 
pattern with its maximum over the area of mountain crest where loamy clay and 
peat can be found. The corresponding development of lower pressure in this area 
then led the parcels to converge and to trigger clouds. 

As evident from the near-surface wind field (Figure 7.11b), up-slope winds on the 
slopes of the Harz have developed already at 0830 UTC as well as up-valley winds 
in some of the valleys which extend from the foreland to the inner centre of the 
Harz, like the one in the southwest. 

Figure 7.12 shows the backward trajectories of the parcels with a starting posi- 
tion at cloud base (Æ 1700 m) at 1110 UTC and then integrated backwards until 
0700 UTC. The 3-D view of the backward trajectories of the investigated con- 
vective system shows that the air parcels originate from near-surface positions of 
the surrounding flat forelands of the Harz in the north and south. As transported 
upwards along the slopes of the Harz by upslope winds, they reach the top of the 
Harz, as evident from the position of the trajectories at 0830 UTC. Additionally, 
one branch of trajectories comes along the valley in that is positioned in the the 
southwest slope of the Harz (Figure 7.1 1a) and finally feeds the convective cloud. 
At that time, all trajectories are concentrated within an area with a diameter of 
about 15 km (Figure 7.1 1a). 


7.2. Secondary triggering 


7.2.1. Cold pools 


As the name suggests, cold pools are the evaporatively-cooled regions of down- 
draught air. They typically form in the boundary layer as a result of convective 
precipitating systems. The evaporating rain and hydrometer loading creates the 
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Figure 7.11.: Surface sensible heat flux (colour shaded) and the trajectory points (circles 
with air temperature in the colour shade) (a), wind convergence (colour 
shaded) and wind vectors at 10 m (b) at 0830 UTC, in Aj56,, of case over 
A2, dated 09 June 2018. 


1.4? 


density anomaly (Torri et al., 2015). In this way, the lateral expansion of this cold 
and dry region is governed by the density current dynamics (Charba, 1974). This 
cold air creates a region of enhanced surface pressure which further generates a 
horizontal pressure gradient force that drives the outward flow of cold pools and 
thus forms an outflow boundary called as gust front. Therefore, the ideal para- 
meter for the identification of these density currents (cold pool) is the density or 
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Figure 7.12.: Backward Lagrangian trajectories using LAGRANTO in 415672 of case over 
A2, dated 09 June 2018. The black crosses denote the starting points at the 
cloud base at 1010 UTC and the red circles denote the end points of traject- 
ory integration until 0700 UTC. The colours along the trajectories indicate 
the air temperature (K) and the colour in the xy plane denotes orography. 


virtual potential temperature (Emanuel, 1994). This is also a prognostic variable 
used in ICON and a model output field, which is defined as: 


R, 
Nar Y æ (7.7 


6, = O[1+( 
Ra k#£v,d 


where 6, is the virtual potential temperature, 0 is potential temperature, Rg and R, 
are the gas constant of dry air and water vapour, having values 287.04 Jkg"!K! 
and 461.5 J kg"! K“! , respectively. q, and g, are the specific humidities of water 
vapour and condensate water species, respectively. Cold pools play a crucial role 
in the convective organisation of clouds and therefore, modulating the coupling 
between the CBL and free troposphere. 

Moreover, cold pools aid the formation of deeper, wider and more buoyant clouds 
in association with higher precipitation rates, which further lead to the creation 
of stronger cold pools (Bóing et al., 2012). According to Weisman and Rotunno 
(2004), cold pools also play a crucial role in the maintenance of long-lived squall 
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Figure 7.13.: The areal mean SMRC for SMI in black, TAI in green, orography (Orog) in 
blue, net radiation (NetRad) in pink, and horizontal wind at 10 m (Wind10m) 
in cyan with respect to H (a) and E (b), respectively for A625, of case over 


A2, dated 09 June 2018. The vertical dotted line shows the onset time of 
first cloud in this investigation region. 


lines. They presented the analysis of the simplified two-dimensional vorticity - 
stream function model, which confirms the role of shear in the modulation of flow 
around the cold pools. This study proposed that the buoyancy gradient at the cold 
pool edge generates the horizontal vorticity at the gust front, which can easily 
interact with the vorticity generated by the environmental wind shear and result 
into vertical updraughts. (Rotunno et al., 1988; Weisman and Rotunno, 2004). 
This sort of lifting is called as mechanical forcing, where the near-surface ambient 
air is lifted by the dense cold pool air mass, and thus initiates new convection. 
The second mechanism through which cold pools can initiation new convection 
is thermodynamic forcing. Necessarily, cold pools can trigger convection also in 
an environment of low vertical wind shear. The decaying cold pool could leave 
moisture behind via its precipitating downdraught, and this accumulated moisture 
may also provide enough buoyancy so that a small forcing could also trigger a new 
convection (Tompkins, 2001). This process is termed as thermodynamic forcing. 
Thermodynamic forcing creates the region of favourable convection initiation, and 
the mechanical forcing is needed in order to lift the parcels so that it can overcome 
the CIN and reach the LFC (Torri et al., 2015; Jeevanjee and Romps, 2015). 

Though the cold pools are itself the structures developed after the precipitating 
system, cold pools can also trigger new convection through other different kinds of 
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forcings, for example, by intersecting cold pools itself. Therefore, cold pools res- 
ult in "primary-generation" precipitation initiating "secondary-generation" shal- 
low to deep clouds or resulting in the organisation of clouds. The convective ini- 
tiation can occur preferably at the region of localised differences in the turbulent 
fluxes, low-level diabatic heating and moisture transport. This develops a pressure 
gradient over the area and therefore locally enhanced low-level wind convergence. 
According to Atkins et al. (1998) and Bluestein et al. (1990), the local deepening 
of a moist layer tends to occur at the intersection point of a dryline and an outflow 
boundary (cold pool). The outflow carries the gust front at its leading edge which 
lifts the air parcels ahead of it to the LCL and acts as a primer element for the 
formation of deep convection. 


Main characteristics of cold pools 


As discussed in Section 6.3, after self-aggregation of shallow cumulus clouds, 
precipitation occurred in A156, in the east of Berlin and eventually resulted in cold 
pool formation. Following the general characteristics of cold pools as described 
in Torri et al. (2015); Drager and van den Heever (2017), Figure 7.14 shows the 
spatial pattern of 0, anomaly, q, anomaly, and vertical velocity at Æ1 km, giving 
an example of the cold pool formation in 41567౫ for case simulated over Al, dated 
26 July 2012. The characteristics are: 


- The innermost region is colder (by IK), negatively buoyant and drier than 
the average (or boundary layer air) (Figure 7.14a,b). 


- Cold pool edges characteristically pose negative temperature anomaly and 
enhanced water vapour than the average (Figure 7.14a,b). 


- Additionally, the cold pools are associated with gusty winds (Figure 7.14c), 
so that the horizontal outflow bounded with the gust front would propagate 
with a speed which would be proportional to the density difference, therefore 
buoyancy, between the cold pool and the environment air. 


Superposition of orographically-induced wind and cold-pool outflow 


In this section, an example of convective precipitation generated through low- 
level wind convergence is discussed for Asm of case A2, dated 09 June 2018. 
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Figure 7.14.: Spatial pattern of virtual potential temperature anomaly 6,’ at 10 m (a), spe- 
cific humidity anomaly q,’ at 10 m, and vertical velocity at ~1 km (c) at 
1500 UTC in A156, of case over Al, dated 26 July 2012. 
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Figure 7.15.: Orography (a), cloud cover at 1250 UTC (b) and the resulting accumulated 
precipitation at 1500 UTC (c) in A156m of case over A2, dated 09 June 2018. 


The convection developed in the southwest of Harz mountains (Figure 7.15). This 
convective precipitation 1s in continuation with the above discussed up-slope wind 
example (Section 7.1.2). The diurnal cycle of the accumulated mean precipitation 
of this day shows two periods with precipitation events (Figure B.10). The precip- 
itation during the first period (1100 to 1400 UTC) is the result of clouds formed 
over the Harz mountains due to primary triggering processes (Section 7.1.2). This 
is followed by precipitation in the second period (1500 to 1800 UTC). It is the 
result of clouds triggered in the southwest of Harz due to the superposition of 
primary and secondary triggering mechanisms. This region comprises of mod- 
erately elevated areas with some valleys (Figure 7.15a). The clouds are formed 
at 1200 UTC (Figure 7.15b). With ambient east-southeasterly wind, these clouds 
precipitated out around 1320 UTC near the southeastern tips of the valleys of this 
region (Figure 7.15c) and resulted in precipitating cells with a maximum precipit- 
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ation core of 2:22 mm. To understand the triggering mechanisms of these clouds, 
the low-level wind field and wind convergence are shown in Figure 7.16. There 
are three branches of wind causing wind converging in the southwest of the Harz 
mountain. The first branch is the flow deflected by the Harz mountain itself, com- 
ing from southeast. To determine why the flow on this day goes around instead of 
over the mountain can be explained by the Froude number Fr = U /Nho (Smith, 


1989), where U is the horizontal wind speed, ho the height of the obstacle and the 
g 90, 

8, dz 
the flow has enough kinetic energy to surmount the mountain, while Fr < 1 shows 


Brunt- Vaisala frequency N = A Froude number of Fr > 1 indicates that 
blocking of the flow so that the flow goes around the mountain. The mean Fr for 
the layer (~240 m a.g.l.) at 1000 UTC is Æ 0.41. This explains why the flow goes 
around the Harz in the lowest layers rather over the mountains, and this deflec- 
ted southeasterly wind forms the first branch of the low-level wind convergence 
(Figure 7.16). Most likely, the deep convective cell over the Harz also supports 
a deflection of the approaching southeasterly wind. The second branch consists 
of the valley winds generated in the valleys located in the west of Harz moun- 
tains. These valley winds are further assisted by the wind flowing around the Harz 
in the north and turning into a northwesterly wind at the western tip of the Harz 
(Figure 7.16). Lastly, the third branch generating the wind convergence is the 
cold pool outflow. The deep convective cells which formed over Harz mountains 
subsequently resulted into cold-pool formation. The gust front associated with 
the cold pool assisted the wind convergence in the region southwest of the Harz 
mountain, which ultimately triggers the cloud formations. The contribution of the 
three flows contributing to the wind convergence can also be seen evidently in the 
backward trajectories calculated for these clouds. The backward trajectories are 
initiated from cloud base at 1250 UTC and integrated backwards until 0900 UTC 
(Figure 7.17). The trajectory lines follow the southeasterly wind flowing around 
the Harz, the northerly valley winds, and the cold air parcels originated from the 
ridges of Harz mountain. 

The three branches can be further investigated through the positioning of these 
trajectory points relative to the surface sensible heat distribution (Figure 7.18). 
At 0900 UTC the first set of warm air parcels along the southern flank of Harz 
can be seen which constitutes the deflected east-southeasterly wind in this region 
(Figures 7.16 and 7.18a). The second set is the trajectory points located along the 
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Figure 7.16.: Wind convergence (colour shaded) and wind vectors at ~200 m at 1230 
UTC, in A156m of case over A2, dated 09 June 2018. 
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Figure 7.17.: Backward Lagrangian trajectories using LAGRANTO in 15672 of case over 
A2, dated 09 June 2018. The black crosses denote the starting points at the 
cloud base at 1250 UTC and the red circles denote the end points of traject- 
ory integration until 0900 UTC. The colours along the trajectories indicate 
the air temperature (K) and the colour in the xy plane denotes orography. 
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northwest - southeast oriented valley axis and which are colder and denser, travel- 
ling in the southeast direction (Figures 7.16 and 7.18a). As a result of precipitation 
occurrence over the ridge of Harz, the cold pool forms and at 1200 UTC the cold 
and moist parcels originated from the cold pool itself assisted the convergence line 
in the southwest of Harz (Figures 7.18a and 7.18b). The low-level wind conver- 
gence zones originate from three different sources attributed to the triggering of 
these clouds. The wind flow over this region 15 summarised in Figure 7.19. At 
0900 UTC the ambient east-southeasterly and the northwesterly up-valley wind 
already exist with an average wind speed of 3 and 4 ms !, respectively. At 1200 
UTC, the strong northeasterly wind (>8 ms!) originated from the cold-pool out- 
flow, the northwesterly wind which indicates the up-valley wind, and the ambient 
east-southeasterly wind exist. This superposition of lee-side wind convergence 
and cold-pool outflow in the southwest of Harz resulted in cloud formation over 
this region. These thermally-induced and orography-based triggering mechan- 
isms over A2 and A3 are identified in all of the LES grid spacings. But they show 
slight differences in the location of triggering of clouds depending upon the un- 
derlying land-surface resolution. As such, the dominating triggering factors over 
A2 like the soil type inhomogeneity, the valleys are almost resolved at 5000 m and 
1250 m. Therefore, the modification in the land-surface resolution over this scale 
has a small effect on the triggering mechanisms (not shown) and the precipitation 
patterns (Figure B.14). In particular, the coarsening of land-surface resolution 
from 156 m to 1250 m results in a relative precipitation difference of 10 % and 
up to 25 % with 5000 m (Figure 6.9b). However, the resulting precipitation var- 
ies up to 75 % across the given LES grid spacings (Figure 6.9a) (lower for A156; 
and higher for A1250m) as the clouds are smaller in size and therefore, exposed to 
stronger evaporative cooling (Section 6.3). 


Lesson drawn 


The above-mentioned thermally-induced and orography-based triggering mechan- 
isms over different areas are identified in all of the LES grid spacings. Differences 
depend on whether the relevant land-surface heterogeneities are resolved by the 
model grid spacing or not. Therefore, the land-surface resolution determines if 
precipitation is triggered and how the precipitation patterns look like. One very 
process of precipitation initiation is associated with the superposition of different 
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Figure 7.18.: Surface sensible heat flux (colour shaded) and the trajectory points (circles 
with air temperature in the colour shade) at 0900 UTC (a) and 1200 UTC 
(b) in Aj56 Of case over A2, dated 09 June 2018. 


trigger mechanisms, as described in the example for the Harz mountain, where 
orographically-generated convergence and a cold-pool outflow are responsible for 
convective precipitation. 


101 


7 Identified triggering mechanisms 


Ny 


3 


wind speed (m s^) 
EE [00:1.0 అః [2.0:3.0 mm [4.0:5.0) త [6.0:7.0 Mm [8.0 : 9.0) 
m [10:2.0 mem [3.0:40) ఆఅ [5.0:6.0 ఆఆ [7.0:8.0 mm [9.0 : inf) 
Figure 7.19.: Wind rose diagram for wind speed and direction over a region of conver- 
gence zone in the south west of Harz in 1567౫ of case over A2, dated 09 
June 2018. 
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ICON provides an advanced tool to study its modelling skills concerning convect- 
ive precipitation in the LEM mode. This investigation is done within the frame- 
work of the HD(CP)^ project. More specifically, as part of the synthesis module 
(S4: land surface heterogeneity) of HD(CP)? . The focus of this study is to invest- 
igate the impact of model grid spacing and land-surface resolution on clouds and 
precipitation using ICON-LEM simulations and to determine the processes which 
result into the differences across model grid spacing (A,). It furthermore explores 
the capability of ICON - in dependence of grid spacing and land-surface resolution 
- to simulate those processes, which are responsible for the triggering of convect- 
ive precipitations. At first, the geographical areas with a hierarchy of orography 
and associated deep convection Germany-wide are selected. These three regions 
are 1) the flat terrain near Berlin (Al), 11) the isolated mountain range in the central 
part of Germany, called as Harz mountains (A2), and 111) the complex terrain, the 
Black Forest mountains (A3). Six suitable days with low synoptic forcing, two 
for each area, are selected. After selecting the areas and suitable cases, the ICON- 
LEM simulations are designed with six model grid spacings in a nested setup: 
As000m; A2500m; A1250m; A625 m» A312m and Aisn with the same model configura- 
tion. A5000m and A2soom use NWP turbulence and diagnostic cloud cover schemes, 
whereas A1250m to A156m use LES turbulence and grid-scale cloud cover schemes. 
Those simulations along the line from Asgo9 to 615677 Where the model grid spa- 
cing and land-surface resolution changes simultaneously are designated as control 
runs. The control run 41567 is named reference run, and all results are compared 
with respect to this simulation. Moreover, the land-surface sensitivity experiments 
at a scale of 1250 m and 5000 m (L_< res >, O « res > and LO « res >) have 
been performed in order to understand the relative impact of model grid spacing 
(An) and land-surface resolution on convective precipitation. That gives in total 
six sets of control runs and 30 sets of land-surface sensitivity experiments. 
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In the first step, the variability of the land-surface heterogeneity of the areas 
across the given grid spacings is quantitatively described. For that, the areas have 
been classified using the area-mean Coefficient of Variation (CV) of orography, 
soil type and TAI resolved across the selected range of model grid spacings for all 
areas (Al, A2 and A3). The results show that TAI values have the strongest vari- 
ability than orography and soil type. That means TAI could pertain a considerable 
impact on the spatio-temporal behaviour of the locally-induced precipitation. 
Concerning the first science question, 1.6. the impact of model grid spacing and 
land-surface resolution on clouds and precipitation, the results are based on area- 
mean values calculated for the three geographical regions with different degree of 
orographic complexity. The results show that the areal average of accumulated 
precipitation for most of the cases decreases systematically across the LES grid 
spacings from Aj259m to Ajs6m. More specifically, the differences of the areal 
mean precipitation between the control runs (A1250m, A625m, A312m) and the refer- 
ence run normalised by the reference run can be quite considerable, 1.e. the values 
are in the range of -26 to 400 % with the 75th percentile of 155 %. It is also found 
that the precipitating cells in Aj259,, are typically intenser in comparison to the 
other finer LES model grid spacings. Additionally, the onset time of precipita- 
tion can differ by 1 to 2 hours (normally precipitation starts earlier in A1250m than 
A^|565,). This is a consistent model behaviour found in the precipitation patterns 
simulated for all cases. 

The differences in the areal mean precipitation due to the modifications of respect- 
ive land-surface resolutions are considerably small. The relative differences range 
from about -17 to 37 % with the 75th percentile of 7 % with the land-surface res- 
olution of 1250 m and increases to a range of -17 to 49 % and the 75th percentile 
of 22 % with the land-surface resolution of 5000 m. These ranges hold for all type 
of sensitivity experiments (L_< res >, O_< res > and LO_< res >). Comparison 
of the sensitivity runs with reference run also shows that the onset of convection 
nearly occurs at the same time, the precipitation patterns are in the same regions, 
and the intensity of the precipitation cells is very similar. 

The comparison of precipitation between NWP runs (^2500m and Aso00m) with the 
other LES control runs (A156,; to A1250m) shows that both considerably higher (Al 
and A2) or lower (A3) values can occur. For one case in A3, A2500m and Asooom 
runs do not generate precipitation, although precipitation exists in the LES runs. 
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To investigate whether this is due to model grid spacing or turbulence scheme, 
both NWP and LES turbulence closures were applied in an additional 412507 
run. The comparison shows that already the transition from LES to NWP tur- 
bulence parameterization leads to considerably different precipitation behaviour 
and amount. That means the differences between the step from LES (4125077) to 
NWP (Ao»590,,) simulations is caused by both the parameterization scheme and 
model grid spacing. 

The above findings pose the question of responsible factors and processes, 
causing such large differences in precipitation across the model grid spacings (Aj). 
To investigate and answer this follow-up question (second science question), one 
case from Al was selected, which indicates quite substantial differences in pre- 
cipitation between the LES control runs A156,, and Aj259m. Consistent impacts of 
model grid spacing (A;) on the onset of CBL- convection and cloud evolution are 
found. 

Additional insight in the cloud evolution and precipitation of Aj56,, and A559, 
which show a similar onset time of CBL clouds but the earlier onset of precipita- 
tion in Aj250m (at 1200 UTC) than in A156, (at 1400 UTC), is gained by analysing 
the macro-physical cloud properties in combination with the heat and moisture 
budget. It is found that the model grid spacing directly impacts the macro-physical 
cloud properties. For example, at the beginning of cloud formation, a higher num- 
ber of small and more scattered clouds are produced in Aj56,, than Aj259m. In 
the subsequent hours, these CBL-based clouds grew deeper into the free tropo- 
sphere in both of the grid spacings. However, while the thicker clouds in Aj259 
reached levels of 8 km and started to precipitate already at 1200 UTC, in A156; 
the cloud tops remained much lower (approximately 6 km). According to the 
areal mean budget calculations, in A156m in the layer between 3 to 5.5 km the 
microphysics term indicates a temperature decrease due to evaporative cooling. 
A horizontal cross-section of the microphysics term at 4 km shows that the areal 
average evaporative cooling is mainly generated by strong cooling at the edge and 
shell regions of the clouds. This eventually results in the dissolution of most of 
the clouds, particular because this layer of the atmosphere is still quite dry. At 
midday, the dissolution of clouds even suppresses the precipitation in Aj56,, but 
simultaneously enhances the moisture in the corresponding layers where evapora- 
tion occurs. Later on, cloud aggregation is an important factor with respect to the 
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onset of precipitation: as soon as small clouds aggregate, the evaporative cool- 
ing, which caused the dissolution of clouds by entrainment, reduced considerably 
and the clouds could grow deeper into the atmosphere. Subsequently, the cloud 
clusters result in precipitation, too. 

Next, the spatial distribution of precipitation is analysed. This investigation 
can be seen in conjunction with the third objective of the thesis, 1.6. the triggering 
mechanisms in relation to grid spacing and land-surface resolution. Overall, it is 
found that the precipitating patterns of the control runs (from A156, to A1250m) are 
quite similar. This finding also holds for the precipitation patterns between the ref- 
erence and the corresponding land-surface sensitivity experiments, as long as the 
land-surface resolution is < 1250 m because most of the land-surface patterns are 
still resolved. However, for land-surface resolutions of 5000 m, the precipitation 
patterns could be quite different. 

Some of the precipitation patterns found in the aforementioned simulations and 
the responsible triggering mechanisms are then investigated in more detail. The 
identified triggering mechanisms are classified into two categories, the thermal 
and dynamic triggering mechanisms which can also be decomposed into primary 
and secondary triggering, based on the origin of the respective triggering. In 
the first category, the typical thermal or primary triggering algorithm is due to 
thermally-induced circulations, generated by land-surface heterogeneity or oro- 
graphy. The examples presented show the corresponding thermally-induced trig- 
gering mechanisms, like lake-breeze circulations and UHI, especially active over 
the flat terrain. 

Further, the thermally-induced circulations have been effective together with the 
orographically-induced circulations. The superposition of two of the triggering 
processes seems more effective among the identified example, e.g. superposition 
of upslope-, up-valley winds and soil-type inhomogeneity-based circulations over 
the Harz mountains. This led to intensive precipitation over the mountain ridge, 
which eventually formed the cold pools. This further led to the example of the 
secondary triggering process on this day. 

The other category comprises of typical dynamic and secondary-triggering mech- 
anism, for example, triggering from cold pools. In the given example, the trig- 
gering is most effective when orographically-induced lee-wind convergence su- 
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perimposes the convergence zones of the cold-pool outflow. It is shown that the 
secondary triggering does not occur when one of these preconditions are missing. 
The overall comparison of triggering processes in the control runs together with 
the land-surface sensitivity experiments show that the model grid spacing of ~ 
A1250m and finer led to a similar set of triggering mechanisms. However, convec- 
tion which is triggered under high land-surface resolution does not occur when the 
land-surface resolution is too coarse (e.g. 5000 m) to resolve those land-surface 
inhomogeneities, which are responsible for the generation of thermally-induced 
circulations and the subsequent initiation of convection. Therefore, considerable 
differences can be seen between the triggering processes at land-surface resolu- 
tions of 156 m and 5000 m. Additionally, as proved for A156, and A1250m the 
processes which are responsible for further development of clouds, like evapor- 
ative cooling, could be different and thus be responsible for differences in the 
precipitation patterns. These differences are reciprocated to the range of variabil- 
ity in the areal mean precipitation which is in the order of Æ 22 % in total when 
the land-surface resolution is modified from 156 m to 5000 m. 
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A. Acronyms 


ASTER Advanced Spaceborne Thermal Emission and Reflection Radiometer 
BLIDS Siemens lightning information service 

BMBF Federal Ministry of Education and Research 

CAPE Convective Available Potential Energy 

CBL Convective Boundary Layer 

CI Convection Initiation 

CIN Convective Inhibition 

COPS Convective and Orographically-induced Precipitation Study 
COSMO COnsortium for small-scale MOdelling 

CPMs Convection-Permitting Models 

CSIP Convective Storm Initiation Project 

CV Coefficient of Variation 

DC Deep Convection 

DRM Dynamic Reconstruction Model 

DWD Deutscher Wetterdienst 

ECMWF European Centre for Medium-Range Weather Forecasts 
ECMWF-IFS ECMWF-Integrated Forecasting System 


ERA ECMWF Re-Analysis 
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ESA European Space Agency 
EUCLID European Cooperation for Lightning Detection 


EX TPAR External Parameter for Numerical Weather Prediction and Climate 
Application 


FKZ Forderkennzeichen 

HYSPLIT Hybrid Single Particle Lagrangian Integrated Trajectory model 
HWSD Harmonized World Soil Database 

ICON ICOsahedral Nonhydrostatic 

ICON-EU ICON-Europe 

ICON-LEM ICOsahedral Nonhydrostatic - Large Eddy Model 
IHOP International H2O Project 

LAGRANTO LAGrangian ANalysis TOol 

LAI Leaf Area Index 

LAM Limited Area Mode 

LCL Lifting Condensation Level 

LEM Large Eddy Model 

LES Large Eddy Simulation 

LFC Level of Free Convection 

LHR Latent Heat Release 

LSMs Large-Scale Models 

MetUM UK Met Office Unified Model 


METI Ministry of Economy, Trade, and Industry of Japan 
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MLCAPE Convective Available Potential Energy of mean surface layer parcel 
MLCIN Convective Inhibition of mean surface layer parcel 

MPI-M Max-Planck-Institute for Meteorology 

MTSAT-1R Japanese geostationary meteorological satellite 

NASA National Aeronautics and Space Administration, United States 
NICAM Nonhydrostatic icosahedral atmospheric model 

NWP Numerical Weather Prediction 

PDF Probability Density Function 

RADOLAN Radar Online Adjustment 

RBF Radial Basis Function 

RPD Relative Percentage Difference 

RRTM Rapid Radiation Transfer Model 

SMI Soil Moisture Index 

SMRC Standardized Multiple Regression Coefficient 

TAF Topographic Amplification Factor 

TAI Transpiration Area Index 

TERRA-ML TERRA - Multi Layer 

TKE Turbulent Kinetic Energy 
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(b) Orography over A2 
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(c) Orography over A3 
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(f) Soiltype over A3 
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(1) TAI over A3 


Figure B.1.: Composites of modelled orography (a-c), soil type (d-f) and TAI (g-1) across 
the chosen model grid sapcings (Asooom» Ao500m. A1250m» A625m» À312m 
A156m) Over flat terrain (Al), isolated mountain range (A2) and complex ter- 
rain (A3). 
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Figure B.2.: Horizontal distribution of lightning density (a) and observed (RADOLAN 
RW) daily accumulated precipitation (b) over Al on 17 Sept 2017. 
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Figure B.3.: Same as Figure B.2 but over A2 on 17 Sept 2017. 
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Figure B.4.: Same as Figure B.2 but over A3 on 29 May 2017. 
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Figure B.5.: Total accumulated precipitation at 2000 UTC over Al, case: 17 Sept 2017. 
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Figure B.6.: Total accumulated precipitation at 2000 UTC over A2, case: 09 June 2018. 
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Figure B.7.: Total accumulated precipitation at 2000 UTC over A2, case:17 Sept 2017. 
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Figure B.8.: Total accumulated precipitation at 2000 UTC over A3, case: 29 May 2017. 
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Figure B.10.: Time series of total accumulated precipitation for NWP (Asooom; A2500m); 
dashed lines, and LES model grid spacings(Aj259m, A625m, A312m, A1s6m)s 


solid lines, over A2, case: 09 June 2018. 
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Figure B.11.: Time series of total accumulated precipitation for NWP (Asooom; A2500m); 
dashed lines, and LES model grid spacings(A1250m;> A625 m» A312m> Nis6m); 
solid lines, over A2, case: 17 Sept 2017. 
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Figure B.12.: Time series of total accumulated precipitation for NWP (Asooom; A2500m); 
dashed lines, and LES model grid spacings(À1250m À625m; A312m> A156m). 
solid lines, over A3, case: 12 Aug 2015. 
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Figure B.13.: Time series of total accumulated precipitation for NWP (Asooom; A2500m); 


dashed lines, and LES model grid spacings(Aj259m, A625m, A312m, 156m). 
solid lines, over A3, case: 29 May 2017. 
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Figure B.14.: Total accumulated precipitation from A456," control run (a), Ajs56m- 
LO 1250 m run (b), A156m- L_1250 m run (c), Ajs6m- O_1250 m run (d), 
A156m-LO_5000 m run (e), A156,,-L. 5000 m run (f), and A1565,,-O. 5000 m 
run (g) at 2000 UTC over A2, case: 09 June 2018. 
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Figure B.15.: 
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Figure B.16.: 


128 


08 09 10 11 12 13 14 15 16 17 18 19 20 
Time (UTC) 


Time series of total accumulated precipitation from A156,- control run (L), 
A156m- LO 1250 m run (8), A156m- L 1250 m run (+), A156m- O. 1250 m 
run (0), A156," LO. 5000 m run (ఆ). A156,,-L. 5000 m run (+), and Aj56m- 
O. 5000 m run (o) over A2, case: 09 June 2018. 
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Figure B.17.: Same as Figure B.16 but over A3, case: 12 Aug 2015. 
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Figure B.18.: Vertical profile of total moisture advection (solid) and grid-scale latent heat 
flux convergence (dashed) at 1200 UTC in Aj259 of case over Al, dated 26 
July 2012. The horizontal lines with circle and square markers denote the 
areal mean CBL height (zi) and LCL, respectively. 
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Figure B.19.: Same as figure B.18 but for 615677. 


Ag156m 


—— MICRe, 
—— TURBa. 
— ADV, 
— RAD. 
— SUMg, 
— LCL 
= Z| 


Height (km) 


— MICRav 
—— TURBgy 
—— ADVy 
—— SUMa 
—— LCL 
— 1 


ZG es c4 53 IE Ll 273 4.5 4B 
Moisture budget (g kg! h7!) 


Figure B.20.: Vertical profile of heat (a) and moisture (b) budget terms throughout CBL 
and free troposphere in Ajs6m at 1220 UTC for case over A2, dated 09 June 
2018. The horizontal line in orange and cyan colour shows the LCL and 
mean CBL height (zi), respectively. 
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Figure B.21.: Spatial distribution of temperature microphysics tendency (colour coded) 
(MICRg,) in A1250m (a) and A156,, (b) at ~ 6000 m height a.g.l. at 1220 
UTC for case over A2, dated 09 June 2018. The black contours are the 
isolines for cloud liquid water content of 0.1 gkg" at that level and the 
brown contours encircle the areas of lakes. 
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Figure B.22.: Orography (a), soil type (b), cloud cover at 0920 UTC (c), and the resulting 
accumulated precipitation at 1100 UTC (d) in A625,m of case over A3, dated 
12 Aug 2015. 


131 


B Figures 


wW 
e 
e 


N 
Ne) 
Ne) 


N 

Ne) 

oo 
) 


N 
Ne} 
SI 


N 
Ne) 
[93] 


N 
Ne} 
Ui 


00 


Temperature (K 


N 
Ne} 
> 


00 


N 
Ne} 
w 


8.65 
Longitude (°) 8.75 g.g 48.6 


Figure B.23.: Backward Lagrangian trajectories using LAGRANTO in A625,m of case over 
A3, dated 12 Aug 2015. The black crosses denote the starting points at the 
cloud base at 0920 UTC and the red circles denote the end points of traject- 
ory integration until 0600 UTC. The colours along the trajectories indicate 
the air temperature (K) and the colour in the xy plane denotes orography. 
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Figure B.24.: Surface sensible heat flux (colour shaded) and the trajectory points (circles 
with air temperature in the colour shade) at 0820 UTC, in A625m of case 
over A3, dated 12 Aug 2015. The lakes and water bodies are indicated by 
surface sensible heat flux value (H) of ~-20 W m". 
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Figure B.25.: The areal mean SMRC for SMI in black, TAI in green, orography (Orog) 
in blue, net radiation (NetRad) in pink, and horizontal wind at 10 m 
(Wind10m) in cyan with respect to H (a) and E (b), respectively for A625m 
of case over A3, dated 12 Aug 2015. The vertical dotted line shows the 
onset of first cloud in the investigation region. 
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C. Budget terms 


C.1. Heat budgets 


Heat budgets physically represent the conservation of heat. According to Stull 
(1988), it is based on the First Law of Thermodynamics which basically describes 
the conservation of enthalpy. The different components represent the processes 
contributing to instantaneous heating (or cooling) of the atmosphere. This in- 
cludes the contributions from sensible and latent heat energy transport. The basic 
heat conservation equation 1s: 

90 00 8 1 0905 L,E 


— + U, — 


àr t Ua న PC, D 


SUMe, ADVe, DIFFg, | RADe, _ MICR9, 


where U; is the wind vector having three components (u,v, w), Ve is the thermal 
diffusivity, Cp is the specific heat for moist air at constant pressure, which is re- 
lated with the specific heat of dry air according to the equation C.2, Q5 is the 
component of the net radiation in the jh direction, E is the mass of water vapour 
per unit volume per unit time which is being created by a phase change from li- 
quid to solid phase, and Lp is the latent heat associated with the phase change of 
E. Depending upon different phase changes, the values for latent heat at 0° are 
Ly = 2.50 x 10°J/kg (gas liquid), Ly = 3.34 x 10°J/kg (liquid solid), and 
L, = 2.83 x 10°J/kg of water (gas — solid). 


Cy = Cy4(14- 0.84q) (C.2) 


In equation C.1, SU Me, denotes the local storage of heat. The other terms rep- 
resent the heat advection by mean wind (ADVs,), the mean molecular conduction 
of heat (DIF F¢,) (which is relatively small), and the mean net body source term 
associated with radiation flux convergence (RADg,) respectively. MICRg, repres- 
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C Budget terms 


ents the microphysics contribution mainly the LHR due to condensation of water 
vapour and also evaporation of clouds and rain. Now, in order to represent the 
mean and turbulent parts of the designated variables in the model, the dependent 
variables in the equation C.1 is expanded into mean and turbulent parts. After ap- 
plying the Reynolds averaging and putting the turbulent advection term into flux 
form gives equation C.3. 


006 | ..90 8 1 005 LE  Quj'0" 

zd m =U "MOL SS sh E= = 

dt ð j 9x; PC, Ox; PC, Ox; (C.3) 
= \  —u— - M 
SU Me, ADV, DIF Fo, RADo, MICRe, 70800, 


The equation C.3 is similar to the equation of basic conservation of heat C.1, if 
we replace the instantaneous heat with the mean forcing in the following equa- 
tion, and then add the turbulence term (TURBg,) at the end. This additional term 
represents the convergence of turbulent heat flux (unresolved). 


C.2. Moisture budgets 


The formulation of moisture (q,) budget terms are shown here as per Stull (1988). 
Following the same steps as C.1, the conservation of moisture can be written (as- 
suming incompressibility) as, 

Ody ay Sw , E 


dqy 
p e = C.4 
Ot Ox; Vav ox; B Pair P Pair ( ) 


SUM, ADV, DIFF, MIS, MICR,, 


The terms SU M,, represents the net water vapour storage tendency, ADV}, denotes 
water vapour advection, DIF F}, is the diffusion terms (which is again very small 
relative to the other terms). MIS, represents the convergence of falling liquid or 
solid water, for example precipitation, which is not advecting with the wind. It is 
a net body source term. Term MICR,, represents the microphysics contribution 
due to the conversion of solid or liquid into water vapour and vice versa, basically 
microphysical water vapour source (or sink) rates. Now, splitting the dependent 
variables into mean and turbulent parts, averaging the equation and apply Reyn- 
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old’s averaging rules, and using the turbulent continuity equation to transform the 
turbulent advection term into flux form, the resultant equation would be: 


9% — OG, 9% I Qu jq," 
zip o Ran త ue = (C.5) 
Ot Ox; OX; Pair Ox; 
~ — € / 
SUM,, ADV,, DIFF, MICR,  TURB,, 


Similar to before, equation C.5 represent the respective terms except for the addi- 
tion of the turbulence term at the end (TURB), which represents the convergence 
of turbulent moisture flux. In the atmosphere, the turbulent diffusion is on sev- 
eral orders of magnitude more than the molecular diffusion. In this study, the 
area-averaged vertical profiles have been analysed to understand the exchange of 
heat and moisture between surface and boundary layer, and further the exchange 
between CBL and free troposphere. The terms in equation C.3 and C.5 are calcu- 
lated in the model using the given prognostic variables. For more theoretical and 
formulation details see Stull (1988); Langhans et al. (2012). 


Grid-scale and subgrid scale flux 


Another important formulation used in this study involves the decomposition of 
the total advection (ADVa,) into a horizontal (HADVg,) and vertical (ZADV6, ) 
parts. Starting with the equation C.5: 


90, 900, 08, 9 0, 


= C.6 
a er (C.6) 
— — 7 
ADVo, HADVa, ZADVe, 


From the continuity equation which states the conservation of mass, we have: 


OU; _ 


Ds —0 (C.7) 


If we multiply the continuity equation C.7 by € (which denotes any variable, such 
as temperature, specific humidity), we get: 2d = Q0. As an example for 6,: 
J 


a (C.8) 
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Since this term is equal to zero, adding it to the total ADV term of the equation C.6 
will be also equal to zero. 


ADV, = Uj + 032, (C.9) 
90,U; 

ADVa, = — 7 (C.10) 
J 


Equation C.10 is also called as the flux form of the advection term or an expression 
from anelastic approximation. Since any resolved variable can be decomposed 
into horizontally averaged part and a resolved deviation from the spatial average. 
For example 6, can be decomposed as: 


0, = 0,4- 0// (C.11) 


This way, the 0, and vertical wind velocity w in the terms ZADVg, of equation C.6 
can be decomposed as: 


9 (0, + 6) 


00 
ZADVo, =w— = (W+w’) 5 
గ 


RP (C.12) 


After simplifying and applying Reynold’s averaging, which states that |(A.B) = 

AB+a'b' 

owe, n dw’ 0,/ 
Oz dz 

ల. RE 


ZAm_divg ZAGS_divg 


ZADVo, = (C.13) 


where ZA; jive and ZAgs,ivg are the mean vertical- and grid scale- sensible heat 
flux divergence, respectively. The turbulent components of ZADVa, has been used 
extensively in the analysis of turbulent exchange processes. The grid-scale sens- 
ible heat fluxes, ZADVgs have been calculated using the model output fields of w 
and 0,. 


ZADVgs = (w = w)(0, = 0,).p (C.14) 
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! virtual potential temperature 


where w, 6,, and p are vertical velocity in ms” 
in K and density in kgm ^? at each grid points, respectively. w and 0, are the 
areal mean over the considered domain. Further, subgrid sensible heat fluxes are 
calculated using the model output fields, turbulent diffusion coefficients for heat 
Kr (‘tkvh’), 6. and the geometric height of model levels above sea level (NN) 


z (‘HHL’) (Neske, 2017). 


dé, 


ZADVsgs = —Kn 
dz 


(C.15) 


2 


where K; is in m s=! and 0, is in K, and z is in metres. 
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